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ABSTRACT 


This report presents an analytical investigation of the effect of 
unmodeled measurement system errors on the accuracy of aircraft stability 
and control derivatives identified from flight test data. Such error 
sources include biases, scale factor errors, instrument position errors, 
misalignments, and instrument dynamics. Output error identification 
algorithms that tend to minimize quadratic functions of the difference 
between actual and modeled aircraft trajectory measurements are studied. 

Two techniques - ensemble analysis and simulated data analysis - are 
formulated to determine the quantitative variations to the identified 
parameters resulting from the unmodeled instrumentation errors. The 
parameter accuracy that would result from flight tests of the F-4C air- 
craft with typical quality instrumentation is determined using these 
techniques . 

It is shown that unmodeled instrument errors can greatly increase 
the uncertainty in the value of the identified parameters. Some improvement 
can be made to the identification accuracy by treating the error sources 
as unknown parameters and identifying them along with the stability and 
control derivatives. Additional accuracy improvement can be obtained by 
choosing elements of the identification cost algorithm's function weighting 
matrix so that the sensitivity to the dominant error sources is reduced. 

Computation of the sensitivity matrix of aircraft parameter deviations 
to individual instrumentation error sources is made to enable determining 
what statistical variations the identified parameters will have due to each 
of the error sources. This sensitivity matrix is also used to specify 
instrumentation quality necessary for obtaining aircraft parameters to a 
desired level of accuracy. 

General recommendations are made of procedures to be followed to insure 
that the measurement system associated with identifying stabilty' and control 
derivatives from flight test provides sufficient accuracy. 

iii 



CONTENTS 

Page 

I. INTRODUCTION 1 

II. DEVELOPMENT OF ERROR ANALYSIS TECHNIQUES 

2.1 Modified Newton-Raphson Parameter Identification Process 5 

2.2 Linearized Aircraft Equations of Motion and the 

Me as ur ement s 11 

2.3 Effect of Measurement Errors on the Identification Process.... 15 

2.3.1 General Instrument Error Models 16 

2.3.2 Particular Errors Studied 18 

2.4 Ensemble Analysis of Measurement Error Effects 26 

2.5 Simulated Data Analysis of Measurement Error Effects 32 

2.5.1 Linear Analysis 33 

2.5.2 Nonlinear Analysis 36 

2.6 Implementation of the Analysis Techniques 38 

III. STUDY OF THE EFFECTS OF MEASUREMENT ERRORS ON PARAMETER IDENTIFICATION 
ACCURACY 

3.1 Current Flight Instrumentation Accuracy 39 

3.2 Model of the F-4C Aircraft 41 

3.3 Effect of Longitudinal Measurement Errors 43 

3.3.1 Basic Instrumentation Error Effects 43 

3.3.2 Effect of Changed Input and Data Span 50 

3.3.3 Effect of Fewer Instruments 54 

3.3.4 Effect of Changing the Algorithm Weighting Matrix 59 

3.4 Effect of Lateral Measurement Errors..... 63 

IV. SUMMARY AND CONCLUSIONS 74 

V. RECOMMENDATIONS 77 

APPENDIX A - SUMMARY OF EQUATIONS CODED IN THE ANALYSIS PROGRAM 79 

REFERENCES 90 


iv 



FIGURES 


Page 

1. The Output Error Identification Process 6 

2. Effect of Measurement Errors on the Identification Process 19 

3. Model of Sampling and Quantization Effects on 

Measurement Signals 24 

4. Longitudinal Reference Trajectory and Elevator Deflection 

Input 44 

5. Comparison of Standard Deviations of Longitudinal Motion 

of Parameter Estimates for Baseline Instrument Errors 47 

6. Effect of Center-of-Gravity Postion Uncertainty on the 

Estimation Accuracy of the Parameter Z6e 53 

7. Effect of Various Instrument Biases on the Estimation 
Accuracy of the Parameter Mu. Also Shown is the Effect of 

Not Using the Angle-of-Attack Vane Data 58 

8. Effect of Pitot Tube Bias on the Estimation Accuracy of 

the Parameter Xw 60 

9. Lateral Reference Trajectory and Input of Aileron and 

Rudder Deflections 65 

10. Comparison of Standard Deviations of Lateral Motion Parameter 

Estimates for Baseline Instrument Errors 67 


v 



TABLES 


Page 

1. Standard Deviations of Longitudinal Instrument Errors 40 

2. Reference Stability and Control Derivative Parameters for 

the F - 4C Aircraft 42 

3. Standard Deviations of Parameter Estimates Due to Instrument 

Errors: Longitudinal Equations of Motion ■ 46 

4. Longitudinal Parameter Sensitivity to Instrument Errors 

When No Biases Are Estimated 49 

5. Longitudinal Parameter Sensitivity to Instrument Errors 

when Initial Conditions and Biases Except are Estimated 51 

6. Dominant Random Error Sources for Longitudinal Parameter 

Identification 52 

7. Effect of Doubling the Time Span of the Measurement Sequence 

with Constant Input Energy 55 

8. Effect of Instrument Removal on Parameter Deviations for Baseline 

Set of Errors and Longitudinal Equations of Motion 57 

9. Effect of Changed Noise Variances on Total Parameter Estimate 

Deviations 62 

10. Standard Deviations of Lateral Instrument Errors 64 

11. Standard Deviations of Parameter Estimates due to Instrument 

Errors ; Lateral Equations of Motion 66 

12. Lateral Parameter Sensitivity to Instrument Errors when No 

Biases are Estimated 69 

13. Lateral Parameter Sensitivity to Instrument Errors when Initial 

Conditions and Biases Except b, are Estimated 70 

<P 

14. Dominant Random Error Sources for Lateral Parameter 

Identification 71 

15. Effect of Instrument Removal on Parameter Deviations for Baseline 

Set of Errors and Lateral Equations of Motion 72 

vi 



SYMBOLS 



c.g. 

D, D(p) 

D(p) 

E . 
noise 


e 


e . . 

ij 


vector of biases affecting measurement of y 

vector of biases affecting measurement of u 
biases in the longitudinal output measurements 


biases in the lateral output measurements 


matrix which accounts for cross-product of inertia 
term in the lateral equations of motion 

center-of-gravity 

control measurement matrix relating u to y 

A 

model of D containing parameter estimates p 

covariance of parameter estimate errors due to the 
output noise w^ 

measurement errors 
elements of the T matrix 


®M 

e R 

e 0* 


nx* 


nz ■ 


mean measurement errors 

random constant measurement errors 

scale factor errors in the longitudinal output 
measurements 


vii 




e* 

r 


F, F(p) 

F (p) 

F 

c 

F 

m 

G> G(p) 

c(i) 

g 

H, H(p) 

I 

Ixx, Izz 
Ixz 
J 
K 

LB, Lp, Lr, 


P* 


scale factor errors in the lateral output measurements 


system dynamics matrix containing some of the 
unknown parameters p 

model of system dynamics matrix containing parameter 

/\ 

estimates p 

diagonal matrix representing inverse of control 
measurements' time constants 

diagonal matrix representing inverse of output 
instruments' time constants 

control distribution matrix containing some of 
the unknown parameters p 

model of control distribtuion matrix containing 
parameter estimates p 

acceleration due to gravity 

* 

model of H containing parameter estimates p 
identity matrix 

moments of inertia about the aircraft longitudinal 
and vertical axes 

product of inertia of aircraft in the longitudinal 
plane 

performance index used to estimate parameters in 
the identification process 

conversion from radians to degrees 

L6a, L6r roll moments due to lateral velocity, roll rate, 

yaw rate, aileron deflection, and rudder deflection 
perturbations 


viii 



Mq, Mw, Mu, M6e 


pitching moment due to pitch angle, vertical speed, 
longitudinal speed and elevator deflection perturbations 


NB, Np, Nr, N6a, N6r 


n 


n(k) 


n , n , n 

x y* z 


n , n , n 
xc yc zc 


P 

A 

P 


P, q, r 


p 


b 


yaw moments due to lateral velocity, roll rate, yaw 
rate, aileron deflection, and rudder deflection 
perturbations 

number of points collected in the measurement sequence 
quantity representing error due to signal quantization 
longitudinal, lateral, and vertical aircraft accelerations 

corrections to accelerometer readings 

vector of parameters to be estimated 
estimated value of the vector p 
roll, pitch and yaw attitude rates 
unknown biases to be estimated 


P 


P 


unknown stability and control derivatives to be 
estimated 


P 


IC 


unknown state initial conditions to be estimated 


Q 

R 

R 

t 

T 


T 


c 


u 

U, V, w 


quantization level 

covariance matrix of measurement noise w^ 

covariance matrix of w , 

ci 

matrix representing scale factor errors and cross- 
coupling errors in the measurement of y 

diagonal matrix representing scale factor errors 
in the measurement of u 

control input vector 

forward, lateral, and vertical velocity perturbations 
indicated value of u due to scaling and bias errors 



U„ 


m 


V 


w 


W, 


W- w , w , 
0 * q * nx * 

w 

nz 


X 


vcg 


Xw, Xu 


x, x 


X cg’ y cg’ Z cg 


Y6, Y6a, Y6r 


Zw, Zu, Z<5e 
z(k) 


value of Uj due to measurement system lags 
measurement value of u corrupted by noise 

Li 

nominal airspeed 
control measurement noise 

contaminating output measurement noise assumed to 
be white 

standard deviation of noise in the longitudinal 
output measurements 

angle-of-attack vane distance from the aircraft 
center of gravity 

longitudinal force due to vertical and longitudinal 
speed perturbations 

aircraft state vector 

components of the accelerometer position from the 
aircraft c.g. 

side force on aircraft due to lateral velocity, 
aileron deflection, and rudder deflection perturbations 

measurement vector of aircraft state and its derivatives 

indicated value of y^, due to scaling, cross-coupling, 
and bias errors 

value of y^ due to measurement system lags 
true value of y 

vertical force due to vertical speed, longitudinal 
speed, and elevator deflection perturbations 

sampled signal 


x 



a > a 
o 


nominal angle-of-attack 


Y nx * Y nz 


V Y r 


Y* , Y- 


Ap 

Sp 

At I 

At 


s 


6e 

6 


ij 


6p 


e , e 
ax az 


e , e 
cgx cgz 


vx 


* 

P 

a 


n 


T 

* 


nominal sideslip angle 

misalignments of the longitudinal and vertical 
accelerometers about the aircraft lateral axis 

misalignments of the roll and yaw gyros about the 
aircraft lateral axis 

misalignments of the roll and yaw angular accelerometers 
about the aircraft lateral axis 

error in p obtained from an individual collection 
of measurement data 

mean value of Ap 

time step of the numerical integration method 

sample time step 

elevator deflection 
Kroneker delta 

the difference between p and p 
errors in the accelerometer location 

error in e.g. location 

error in ct-vane location 

nominal pitch angle 
nominal yaw angle 

correlation of consecutive sampled terms 
standard deviation of n(k) 

time delay in sampling control input 
nominal roll angle 


o 


xi 



NOTATION 


( ) 

( 

(") 
E{ } 


3 ( ) 

3p 


K ) 

i=l 

A( ) 


time derivative of variable 
sample of variable at i £ instant 
estimated value of a parameter 
expected value of a variable 
matrix transpose 
matrix inverse 

gradient of a constant with respect to the vector p. 
For a constant, this is taken to be a row vector. 

For a vector, this is taken to be a matrix with the 
number of columns equal to the order of p and the 
number of rows equal to the order of the vector. 

summation of points from 1 through n 


perturbation of a quantity about the nominal value 
or trim position 

modification of parameter due to cross-product of 
inertia terns 

measured value of aircraft states and their derivatives 


xii 



I 


INTRODUCTION 


The process of determining stability and control derivatives of an 
aircraft from flight test data is called aircraft parameter identification. 

There are several reasons why this process has developed into a very 
important field of endeavor. These include: 

1. Many instances where the prototype aircraft do not have the same 
characteristics as predicted by their wind tunnel models. The 
cost to the United States government due to out-of-control aircraft 
losses has been subtantial ^ . Major cost and safety considerations 
motivate determining ways of obtaining better knowledge of the air- 
craft parameters; 

2. Requirements for better understanding and calibration of wind 
tunnel testing and its relationship to actual flight vehicle 
performance; 

3. The potential of allowing the deeper understanding of aerodynamic 
phenomena and the relationship to vehicle stability; 

4. Requirements for ground-based simulators which are more accurate 
representations of the aircraft in all flight regimes; 

5. Requirements for superior stability augmentation and adaptive 
flight control systems. 

There are three essential elements in the development of more adequate 
methods for identifying aircraft parameters from flight data: 

1. Improved algorithms and computer programs to identify the derivatives, 
their confidence levels (variances).i and related parameters such as 
sensor errors and wind gusts; 
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2. The determination of proper sequences of flight control inputs 
(surface deflections) which will excite all the aircraft response 
modes from which parameters are to be extracted, and methods of 
displaying this information to the test pilot so that he is aware 
of when a suitable maneuver has been executed; and 

3. Adequate instrumentation (the right kind of sensors with necessary 
accuracy) and recording equipment with which to collect the flight 
data. 

This study is concerned with this last point, namely, the establishment 
of what constitutes instrumentation accuracy to enable the collection 
of flight data which is of adequate quality for identifving the aircraft 
parameters to the accuracy desired. 

In general, flight instrumentation is not specified today for the direct 
intention of identifying stability and control derivatives. Rather, its 
intended purpose is for checking aircraft handling qualities and general 
measures of performance. If instrumentation specification is made, it is 
typically based on what is known to be available. Part of the reason for 
this status is that estimating stability derivatives from flight test data 
has only been a secondary activity of companies building aircraft. If 
a problem arises in the handling qualities, the manufacturer may attempt 
to determine the derivatives responsible for the undesirable characteristic 
as an aid to the best design fix; however, generally no full identification 
program is undertaken. Flight simulators are built using wind tunnel 
estimates of stability derivatives, and only corrections for gross dis- 
crepancies are made. 

There have been two notable exceptions ^ to the lack of attention 
given to specifying instrumentation for the direct purpose of extracting 
stability and control derivatives from flight test data. The Technological 
University at Delft, the Netherlands ^ ^ has developed instrumentation 
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systems with digital data acquisition, precision temperature controlled 
electronics for uniform instrument dynamics, and inertial instruments 
in a temperature controlled housing. However, data of individual instrument 
contributions to identification errors have not been collected, nor have 
the individual error effects on particular stability derivatives been 
determined. 

(3) 

LTV Aerospace Corporation has studied instrument error effects on 
VTOL parameter identification accuracy. The LTV work involved repeated 
simulation of the identification process and included random noise error 
sources. A least-squares identification algorithm was used. The large 
parameter estimate errors which are characteristic of least-squares methods 
in the presence of random measurement noise were avoided by including 
"pre-filters" in the data processing procedure. These analog pre-filters 
were implemented on the aircraft to prevent aliasing in the sampling 
process of digital data acquisition. No individual parameter sensitivities 
to particular error sources were reported in their work, so that instrument 
tradeoff judgements couldn't be made. Rather, one instrument set and its 
accuracy level were defined which met the requirements of a particular 
VTOL testing program. 

The purpose of this present study has been threefold: 

1. The development of techniques, algorithms, and a computer program 
with which to assess the uncertainty due to instrumentation errors 
in the accuracy of the aircraft parameters identified from flight 
test data; 

2. The application of these techniques to examine the variation of 
parameters obtained from typical flight tests with typical instru- 
mentation errors; and 
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3. The determination of the general effects of instrumentation 
quality variations, the type of instruments used, and other 
quantities governing the data collection and identification 
process on the identified parameter accuracy. 

This study is a first step in the overall task of specifying and 
providing adequate flight instrumentation for parameter identification. 
The results determine important factors which must be considered and 
procedures which should be followed to insure the measurement system 
is sufficient. 
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II 


DEVELOPMENT OF ERROR ANALYSIS TECHNIQUES 


Techniques are developed in this section to determine quantitatively 
the parameter variations which would result from using an output error 
identification algorithm in the presence of unmodeled instrument errors. 

It is assumed that the identification algorithm is convergent and that it 
tends to minimize a quadratic function of the difference between actual 
and modeled aircraft trajectory measurements. The modified Newton-Raphson 
identification algorithm is specifically used. It is further assumed that 
a single application of this algorithm can determine the major portion of 
the variation in the identified parameter value due to the instrumentation 
errors . 

2.1 Modified Newton-Raphson Parameter Identification Process 

(4) 

The modified Newton-Raphson algorithm is essentially one of several 
output error identification methods which are used. This basic identification 
process is illustrated in Fig. 1. The algorithm's objective is to choose 

A 

parameters p of a mathematical model of the aircraft so that the difference 
between the output measurements of the model and the actual aircraft are 
minimized. With no measurement errors, external disturbances, or model 
structure inaccuracies, the output errors are minimized when the model 
parameters equal those of the aircraft. Output error identification 
methods have the following characteristics: 

1. They require good initial estimates of the aircraft states and 
the parameters; 

2. They give unbiased estimates in the presence of zero mean white 
measurement noise; 

3. They can be used for identifying the parameters of aircraft with 
both linear and non-linear equations of motion; and 

4. They do not work well in the presence of random disturbances to 
the dynamics (process noise) . 
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FIGURE 1. THE OUTPUT ERROR IDENTIFICATION PROCESS 



In this study, the aircraft equations of motion as perturbed 
from the nominal flight path are assumed to be linear with constant 
coefficients, and of the form 


x = F (p)x + G(p)u ; x(0) = x 

o 


( 2 . 1 ) 


where 


A 

X “ 


aircraft state vector 

A 

u = 


control input vector 

F (p) 

A 

system dynamics matrix containing some 
unknown parameters p 

G(p) 

A 

control distribution matrix containing 
unknown parameters 


of the 
the other 


The identification process identifies the parameters of F and G. 


The output y of this system consists of measurements of the elements 
of x and x. It is modeled as a sampled process bv the equation 


y ± = H(p) X;L + D(p)u ] . + w ± (2.2) 

where H and D are other constant-coefficient matrices also containing 
elements of p. The vector w^ is contaminating noise. The subscript i 
indicates that the output is sampled at time i and processed by a 
digital computer. 


The modeled equations of aircraft motion are of the same order 
and form as Eqs. (2.1) and (2.2), and they are represented by 

x = F(p)x + G(p)u ; x(0) = x q (2.3) 


Here, F(p) and G(p) are formed by using the estimated parameters p. The 
simulated output equation is 


y. = H(p) x + D(p)u 


(2.4) 
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If p equals p and x q = x q , the only difference between y.^ and y 
is due to the measurement noise w^ . In the Newton-Raphson identi- 
fication scheme, it is assumed that is a sequence of zero mean white 
noise vectors with the covariance matrix 


E {w. w. } = R6 . . 
ij iJ 


(2.5) 


Furthermore, it is assumed that the elements of w.^ are independent so 
that R is a diagonal matrix. 

The Newton-Raphson identification technique chooses parameters p 
which minimize the performance index or cost function 

n 

(y i ~ V T R_1 (y i ' y i> (2.6) 

i-1 

where n is the number of points collected in the measurement sequence . 

This is done by iteratively applying the equation 


P k+1 P k 


9 Z J 




-1 * 

T 

8 J 

8 p 


(2.7) 


The first partial of J with respect to p is, from Eq . (2.6) 
n 

3J _ oV / - ,T -1 3y 

j- = -22 <y ± - y ± ) R 


3p 

i=l 

The second partial is 


9 2 J 

9^ 


= 2 


i=l 


3y, T 


3y 


3p 


<\2 A 

i . J -1 8 y i 


X gp 1 r" 1 9 P " (y i " V" R ap 


*The notation 9 ( ) refers to taking the partial with respect to the 

9p 

estimated parameter p. 


( 2 . 8 ) 


(2.9) 


8 



This is often approximated by 


9 2 J = 2 
9P 7 



( 2 . 10 ) 


Equations (2.8) and (2.10) are substituted into Eq. (2.7) to yield 


Pk+1 = \ 



( 2 . 11 ) 


Equation (2.11) is the "modified" Newton-Raphson optimization 
technique. It is applied repeatedly to update p until 
Eq. (2.8) approaches a zero value. 


is 


The variance of the estimated parameter vector due to noise 


{ 5p 4pT } noise ‘ 

9 2 J " 
„ 2 


L 3 p J 


( 2 . 12 ) 


where 


6p 


A 

= p - p 


Equation (2.12) is obtained by assuming that the errors due 
to are small so that J is a quadratic surface in the 
vicinity of £ and p. Then, one can write 



(2.13) 
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where y^ - is w^ and y^ is generated using the correct 
parameter p. Thus, from Eqs. (2.12) and (2.13), 


E jfip 6p T j- = E 


( 

r a 2 T i 

-1 /^ T \ T _ 

r 2 i 

_i j 


3 J 

/ 3J \ 3J 

3 J 


1 

9 

. 3 P 

\ 3p / 3p 

1 

<N 

CU 

<ro 

1 



Because 


[$] 


has no noise dependence, this becomes 


E | dp 6p X [ = 


3p‘ 


-1 


(/ _3J \ T _3J 
A 3 P ) 


3 2 J 




-1 


The inner term is expanded to yield 

■fc)’(i - 


5 fe) T R' 1 (y r 5M 

i»l\3p / 1 1 


1 (yj-yj) 1 R" 1 


J-l J J 3p 

Because the measurement noise is assumed to be white, 

E { (y ± ~ 9 ± ) (y. - ?.) T } = R 5. . 

1 1 J J ’ ij 


The double summation reduces to a single summation, and the 
expectation is replaced by R yielding 


E 




R" 1 R R -1 _^i 
3p 


which, is exactly equal to from Eq. (2.10). By 

substituting this result in Eq. (2.15), the desired relation 
(Eq. (2.12)) is established. 


(2.14) 


(2.15) 


(2.16) 


. (2.17) 
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2.2 


Linearized Aircraft Equations of Motion and the Measurements 


It is assumed that the aircraft begins a maneuver from a 
quasi-steady flight condition with a constant airspeed V, 
angle-of-attack a Q , and pitch angle 0 Q . The roll angle 

<P t yaw angle \p f sideslip angle 0 , and the attitude rates 
p, q, and r are all assumed to be initially zero. The 
equations of motion of small perturbations of the aircraft 

( £ 7 \ 

in the longitudinal plane are ' ’ ' 


Here 


> 1 


Aq 


Aw 

ii 

Au 

O 


c T & 

[a6 


0 

0 

-g sin 6 0 /K 
-g cos 0 q /K 


Aq Aw Au 


1 


0 

0 


A0 


0 

Mq 


Mw 

Mu 


Aq 

+ 

M6e 

V cos 

a 0 /K 

Zw 

Zu 


Aw 


Z6e 

-V sin 

o 

3 

Xw 

Xu 


Au 


0 



(2.18) 


and consists of perturbations in pitch, pitch rate, the normal 
component of relative velocity, and the longitudinal component 
of relative velocity. The control A 6e is the deflection of 
the elevator surface about the trim position. The constant 
K is the conversion from radians to degrees. 


If only the short-period motion is to be studied, the equations 


A0 

Aq 

Aw 


0 

0 

-g sin 9 o /K 


1 

Mq 

V cos 


0 


A0 


0 

Mw 


A9 

+ 

M<$e 

Zw 


Aw 


Z6e 

_ 


_ 


_ 



are used . 


(2.19) 



In Eq. (2.18), the unknown parameters to be identified consist 
of 

T A r 

p = [_Mq Mw Mu M6e Zw Zu Z<5e Xw Xu] 

In Eq. (2.19), this reduces to 

m 

p = [Mq Mw M6e Zw Z6e] 

The lateral equations are in the form 

• I T 

Cx = F x + G u 

where C is a matrix which accounts for the cross-product of 
inertia term Ixz. The state 

x = [a 6 Ap Ar A<J>] 


(2.23) 

(2.24) 


consists of perturbations in the angle-of-sideslip, roll rate, 
yaw rate, and roll angle. Then C has the form 



-1 * -1 ’ 

By letting F - C F and G = C G , Eq. (2.22) can be modified 
to the more standard form 

x = C _1 f'x + C _1 g' u = Fx + Gu 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 
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or in full form, 


r • 

A6 


YB 

sin a 

0 

-cos a 

0 

g cos 0 /tf" 


’ ab' 


’Y6 a 

Y6r 

Ap 

_ 

LB* 

Lp* 

Lr* 

0 


Ap 

+ 

'L6a* 

L6r* 

As 


NB* 

Np* 

Nr* 

0 


Ar 


Nda* 

N6r* 

A<j> 


0 

H 

1 

tan 6 

0 

0 


( 

< 

1 


0 

0 


A6a 

A6r 


(2.25) 


The starred (*) quantities are modified from their normal values 
due to C in Eq. (2.24). The control deflections A6a and A<Sr 
are those of the ailerons and rudder, respectively. In 
Eq. (2.25), the unknown parameters are 


T A [ 1 

p = |YB Y6a Y<5r L8* Lp* Lr* L<Sa* L6r* N6* Np* Nr* NcSa* N<Sr* (2.26) 

The seven instruments which are assumed to be available for 
longitudinal measurements are: 

1. pitch attitude gyro (6) 

2. pitch rate gyro (q) 

3. angle-of-attack vane (a) 

4. longitudinal pitot tube or air speed indicator (u) 

5. longitudinal accelerometer (n^) 

6. normal accelerometer (n ) 

z 

7. pitch angular accelerometer (q) 


For the short period equations, the pitot tube and longitudinal 
accelerometer are omitted. 


The lateral instruments are assumed to be 

1. angle-of-sideslip vane (8) 

2. roll rate gyro (p) 

3. yaw rate gyro (r) 

4. roll attitude gyro ($) 
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5. lateral accelerometer (ry ) 

6. roll angular accelerometer (p) 

7. yaw angular accelerometer (r) 



Making the substitutions and fitting the longitudinal measure- 
ments into the form of Eq. (2.4) yields (for linear accelerations 
measured in g’s) 
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The lateral measurements are 


he ^ 

m 


1 

0 

0 

0 

AP m 

m 


0 

1 

0 

0 

Ar 

m 


0 

0 

1 

0 

Ad) 

m 


0 

0 

0 

1 

An 

ym 


VYB 

gK 

0 

0 

0 

• 

Ap 

m 


LB* 

Lp* 

Lr* 

0 

• 

Ar 

m 


NB* 

Np* 

Nr* 

0 

— 


_ 



— 


AB 

Ap 

Ar 

A<t> 


0 

0 

0 

0 

Y6a 

gK 

L6a* 

N6a* 


0 

0 

0 

0 

Y6r 

gK 

L6r* 

N<5r* 


A6a 

A6r 


(2.29) 


Equations (2.28) and (2.29) assume perfect measurements of 
the aircraft state x and the control input u. 


2.3 Effect of Measurement Errors on the Identification Process. 


Often, no other measurement errors except for the white noise 
indicated earlier are assumed to be present in the flight data 
used for identifying aircraft derivatives. Sometimes biases are 
assumed to affect the measurements and these terms are identified 
along with the equation parameters and state initial conditions. 
However, there are many other types of errors which do affect 
the estimation accuracy as will be seen. In this discussion, the 
emphasis is placed on those error sources whose effect can be 
determined by linear analysis. 
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2.3.1 General Instrument Error Models 

First, consider the measurement of the aircraft state. 

For constant value of these outputs the actual indicated 
readings would be of the form 

y I = V B (2.30) 

where 


T 


The diagonal terms in the T matrix represent scaling errors 
while offdiagonal terms represent cross-coupling errors. The 
vector B represents the bias errors. 

The measurements are also affected by the dynamic character- 
istics of the instruments and the recording equipment. The slowest 

(9) 

instrument/smoothing filter combination encountered has a 
natural frequency of 1 cps which is about a factor of 2 
higher than the aircraft dynamics. Therefore, the important 
aspects of the dynamic errors are the phase lag and amplitude 
attenuation of the instruments at frequencies below their 
natural frequencies. These characteristics can be approximately 
simulated by a first order lag regardless of the order of 
the instrument dynamics. The matrix equation representing this 
is 
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(2.32) 


= -Fy + Fy 
m L ml 


y L (0) = y x (0) 


where 

y = "lagged" measurement 

1j 


F^ = diagonal matrix of elements representing one over 
the instruments’ time constants. 

The addition of the random noise for each instrument 
yields the final measurement equation 


y i 


Li 


+ 


w. 

1 


where y^ is the output measurement vector with all errors 
sampled at time i, and w^ is the random output noise vector with 

E {w. } = 0; E {w.w T } = R6^ . 
i i J ij 

In this study, it is always assumed that the random noise is 
correctly modeled; that is, the covariance matrix R is known 
and is correctly used in the cost function J of Eq. (2.6). 


The other source of measurement error is in the recording 
of the control input u by either surface deflection potentio- 
meter or servo measurements. These control measurements are 
also subject to scale factor errors and biases which can 
be represented by the equation 


u T - T u + B 
I c c 

The measurement of u^ is also subject to dynamic effects 
which are again approximated as first order lags by the 
equation 


(2.33) 


(2.34) 
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"l - - F o U L + F c “I i \<-°1 - “!«» (2- 35 > 

Here, 

u^ = "lagged" control 

F * diagonal matrix of one over the time constants of 
the control measurements. 


The actual recorded control input is sampled and is subject 
to noise. It is represented by the equation 


u 


mi 


u. 


Li 


+ w . 
Cl 


where 


(2.36) 


u ^ ~ control measurement vector with errors sampled at time i 


w . = random control noise vector; E (w . } = 0; E {w .w . }= R 6 . . 

ci ci ci cj c ij 

The overall identification process flow diagram changes from 
that depicted in Fig. 1 to that depicted in Fig. 2. In 
the linear analysis which follows, the control measurement 
noise w is ignored. This noise acts as a random disturbance 
to the system dynamics (process noise) and cannot be analyzed 
with linear methods. 

2.3.2 Particular Errors Studied 


Before proceeding to the analysis, a description is first 

presented of some of the error sources which can be studied 

by the preceding equations. The diagonal elements of T, T , F , 

c m 

and F have been explained. B and B are bias vectors. 

c 
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FIGURE 2. EFFECT OF MEASUREMENT ERRORS ON THE IDENTIFICATION PROCESS 





Some of the specific errors which are included in the off- 
diagonal terms of the T matrix include: 

1) a and 8 boom corrections 

2) accelerometer location corrections 

3) misalignments (accelerometers & gyros) 

( 8 ) 

A simplified a boom correction equation is ' 

X 

a = a - vcg q 
V 

where V is aircraft total velocity and X^ is the angle-of- 

attack vane distance from the aircraft center-of-gravity (c.g.). 

If both vane location and c.g. location are precisely known and 

accounted for, there is no error. However, if the actual value 

of X is different from that used in the correction, or if no 
vcg 

correction is made, an error in the measurement results. The 
error in X is thus divided into two parts, the error in vane 
location (e ) and the error in c.g. location (e ). The 
separation of the contributions is made because vane location 
uncertainty only affects the a correction, while c.g. location 
uncertainties affect accelerometer corrections as well. If all 
seven longitudinal instruments are being used as in Eq. (2.28), 
introducing the error Eq. (2.37) into the T matrix Eq. (2.30) 
yields : 


„ e , e 

e 00 = - vx + xcg 

V 

Similar capability is provided for the 8 vane correction errors. 
Other errors which can affect a and 8 readings are due to 
upwash and boom bending . 


(2.37) 


(2.38) 
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Linear accelerometer corrections are necessary when these 

instruments are not mounted at the aircraft center-of-gravity . 

If fx v z 1 are the components of the accelerometer position 
L eg ■'eg cg J 

from the c.g. in aircraft fixed coordinates, then the corrections 
should be ^ ^ 

n xc = (r 2 + q 2 )x cg + (r - Pq)y cg - (q + P r ) z cg (2 ‘ 39) 

n yc = -(r + pq)x cg + (p 2 + r2 )y cg + (P ~ ^ r > z cg 
n zc = (q - rp)x cg - (p + qr)y cg + (p 2 + q 2 ) Z( , g 


These equations can be decoupled into lateral and longitudinal 
parts. If the nonlinear terms are neglected, (valid for p, q, r, 
(expressed in radians/second) which are << 1) , the corrections to 
the longitudinal instruments are 


v 

n = 
zc 


q x 


eg 


(2.40) 


If the value of x and z are in error because of the 
eg eg 

uncertainty in the c.g. position or the c.g. offset of the 
accelerometers is neglected, then the error terms 


57 

(e +e ) 

az cgz 

/Kg 

(2.41) 

67 

(e +e ) 

ax cgx 

/Kg 



appear in the T matrix. In Eq. (2.41), the term 
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errors in the accelerometer location when 
a correction is made. 

distance from c.g. to the accelerometer 
when a correction is not made. 

errors in the knowledge of the c.g. location 

Similarly, the lateral accelerometer has the two errors 


e 56 


- (e + e ) /Kg 

az cgz 6 

(2.42) 

e 57 

= 

(e + e ) /Kg 

ax cgx 


Other elements 
of the gyros and 
the terms 

in the T matrix are due to mounting misalignments 
accelerometers. In the longitudinal equations. 


e 56 

= 

- y /K 
nx 

(2.43) 

e 65 

= 

Y / K 
nz 


appear, where y and y are the small misalignment angles, 
nx nz 

In the lateral equations, the T matrix contains the terms 


e 23 

= 

- V K 

(2.44) 

e 32 

= 

Y /K 
r 


e 67 

= 

-Y-/K 


e 76 

= 

Y-/K 

r 
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which are similarly defined. 


Effects of the above mentioned off-diagonal terms of the T 
matrices of the longitudinal and lateral measurements are 
presented in Section 3. Other errors which could be contained 
in T include angular accelerometer sensitivity to linear 
acceleration and rate gyro mass unbalance. 

Another error source is introduced into the measurements shown 
in Fig. 2 due to the sampling and quantization. This is 
illustrated in Fig. 3a. The errors introduced by this process 
can be duplicated by the addition of a noise source to the 
sampled signal as illustrated in Fig. 3 b. Given the quantization 
level Q and the statistics of the sampled signal z (k) , Widrow 
has developed expressions for the statistics of n(k). 

For all but very course quantization, the distribution of n(k) 
is uniform between -Q/2 and Q/2, and 

E {n(k) } =0 E (n 2 (k)} = Q 2 /12 = a 2 

The error in this approximation is computed based on the relative 
magnitudes of Q and the standard deviation of z(k), (= o z ) , where 
z(k) is Gaussian. When Q > a z ( an extremely course quantization 
level for any airplane measurement system), the error in assuming 
that 


E {n 2 (k) } = Q 2 /12 


is 


2.6 x 10 


-10 



which is very small . 


(2.45) 


(2.46) 


e 
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Perhaps of more interest is the whiteness of the noise sequence 
n(k) , (i.e. , is E {n(k) n(k+l)} = 0?). Widrow also gave an 

expression for this quantity which is 


E {n(k) n(k+l)} = 


2 -(l-p)4iT 2 a 2 /Q 2 

a e z 

n 


where 


E {z(k) z(k+l)} 

P = 2 

a 

z 

For frequent sampling, z(k) will be highly correlated with 

r\j 

z(k+l), i.e. p - 1 so that 1 - p is small (<< 1). 

However, most aircraft measurement systems will have a fine 
quantization level, where cr^/Q is large (>> 1). The net result 
is that it is not clear whether n is white or not. 

As an example, assume that typical numbers for these quantities 

are 

p •" 

a /Q =10 
z 

These yield 


— (1-p ) 4tt 2 a /Q 2 ^ -40 

e z = e 


(2.47) 


(2.48) 
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which is essentially zero. The assumption that the quantization 
adds white noise to the sampled measurements seems reasonable. 

Thus, for the preceding example no special modeling procedure needs 
to be added to include the effect of quantization. 

In summary, there are three types of errors which affect 
the accuracy of instrument measurement of the aircraft motion 
- random noise, random constants such as biases and scale 
factor errors, and mean errors. Mean errors are those terms 
which are known to produce measurement errors but are neglected 
because they are assumed to have negligible effects, (like 
accelerometer offset from the c.g.) Instrument lags can be 
thought of as mean errors with random variations about the 
nominal mean value . 


2.4 Ensemble Analysis of Measurement Error Effects 


As explained in Section 2.1, the modified Newton-Raphson 

identification scheme minimizes the cost function J of Eq.(2.6) by 

repeated application of Eq.(2.11). Convergence on the minimum is 
9 J ~ 

achieved when =0. In this study, it is assumed that the 

true value of p is known. It is further assumed that the instrument- 
ation errors cause the minimum point on the cost function surface 


J to shift a small amount from the true p. If the small error 
assumption is correct, only one application of Eq. (2.6) (with p 
set to p) can determine the shift due to the measurement error on 
the estimate of p. This is the key assumption of the linear analysis 
which is used in this study. The resulting perturbation to the parameter 
vector is: 


<5p 


P " P = 




3J 


3p 


(2.49) 
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where 


p = perturbed parameter estimate due to measurement 
errors 

p = true value of the parameter. 


From Eq. (2.11), this can be written as 


n 

v 

i«yi T 

r- 1 ^4 

-1 

n 3y . T -1 

v 1 J x R 

Zj 

1 3P 

3p j 


Z 9 p 

L i=l 


/ 


i=l 


<yr*i> 


The are the sampled output measurements taken from the aircraft 
(Eqs. (2.1), (2.2) (2.30) - (2.33)) and the y. are the simulated 
output values (Eqs. (2.3), (2.4)) obtained using the measured 
control input. The sensitivity term J i is computed by the 
identification algorithm about the la^lst estimates of p. 

Again, for the linear error analysis, this is the correct value of p. 

3v 

To compute y i in Eq. (2.50) requires integration of 
3p 

k = Fx + Gu ; x(0) = * 

o 

This assumes perfect measurement of the control input u. 

Also, the sensitivities of the states to parameter changes are 
found by integrating 



(2.50) 


(2.51) 


(2.52) 
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where is the parameter vector containing the elements 
presented in Eqs . (2.19), (2.20) or (2.26). If state 
initial conditions are also estimated, the identification 
process integrates 


-d 

dt 




(0) - i 


Then, from Eq. (2.4), the output sensitivity matrix for the 
parameters p is 





/v 

X 


+ 



u 


For the initial conditions, this becomes 


3 £ 

3 Pic 



If output measurement biases are also estimated, the sensitivities 


2 £ _ 

3p b 


i 


must also be included. The total sensitivity used in Eq. (2.50) 
is then 


9y. . 

i A 

w = 





(2.53) 


(2.54) 


(2.55) 


(2.56) 


(2.57) 
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Output measurement errors affect the. value of y_. in Eq. 
(2.50). From Eqs, (2.30) - (2.33), the measurement equations 
can be written as 


y = - F y + F (T{Hx + DuJ + B) 
L mum 


By neglecting the measurement noise temporarily, the output 

y . is the sampled value of y. . The sensitivity of the 
1 L 

error <5p in Eq. (2.50) to an error source e is 


n 





de 


This requires computation of the sensitivity matrix 



The sensitivity of y^^ to a bias element of B is 
approximately 



9 e 3.e 


For an element in T, this is 


9 e 


rs o' 

(Hx + Du) 
d e 


For an unestimated initial condition treated as an error source, 
the sensitivity of y_^ is 


9.y ± 9x 
= H 1 


3;e 


3e 


(2.58) 


(2.59) 


(2.60) 


(2.61) 


(2.62) 
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where — 
3 e 


comes from integrating 


d_ 

dt 




3x 


(0) = I 


(2.63) 


lor an unknown time constant in the matrix F , the sensitivity 

m J 

must be determined by integrating 





+ (Hx + Du) 


7 

' m 3 e 



( 0 ) = 0 


(2.64) 


The results of Eqs. (2.60) - (2.64) are combined into a general 

vector for each error e which affects the output 

3e. 

measurements . 

The sensitivity of parameter estimates due to control input 
measurements errors is of the form 




n 


I 





(2.65) 


This requires knowing the sensitivity of the simulated output 
to control measurement errors e. 
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For control measurement biases, the sensitivity of the control 
input is 


3u 3B 

mi c 

3e 3e 


For scale factor errors, this is 


3u . 3T 
mi _ c 


The effect on the states is found by integrating 


d_ 

dt 



3 x 
3 



from Eqs . (2.3). The resulting effect on the simulated output 

y ± is 

r\ A A rs 

dy. ox. du . 

i _ „ 1 . ™ mi 


Evaluation of the sensitivity to control measurement lags 
requires integrating 



( 2 . 66 ) 

(2.67) 

( 2 . 68 ) 

(2.69) 

(2.70) 
(2.71) 
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The resulting output sensitivity is again found from Eq. (2.69). 

For random constant measurement errors, e , the total 

R 

covariance of the individual parameters being identified is 


E I«P«P T ! total - E |Mp T ! nolse + 



(2.72) 


T 

where E |e^e^| is the covariance matrix of the random measurement 

errors not including measurement noise. The sensitivities 
a (6p) 

9 e come from Eqs. (2.59) and (2.65). The covariance due 
to noise comes from Eq. (2.12). For mean errors, e M , the 
expected error in the parameter is 


E |<5pf 


3(<5p) 

9e M 




(2.73) 


The above error analysis is referred to here as the ensemble 
error analysis. It is valid for small errors which affect 
the measurements linearly. 

2.5. Simulated Data Analysis of Measurement Error Effects 

Sometimes it is useful to determine the effects of instrument 
errors by actually simulating the identification process and the 
measurement data contaminated by errors. If the errors are large 
or if nonlinear errors are to be studied, the one step assumptions 
and linearization which were used in the ensemble analysis method 
may not be valid. Therefore, the simulated data analysis method 
complements the ensemble analysis. 
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2.5.1. Linear Analysis 


There are two common ways in which error analysis is conducted 
by simulation. The first simulates the effect of each random error 
individually, assuming they are independent, and uses the error's 
standard deviation for the magnitude used in the simulation. The 
results of simulating each of the random errors are root-sum-squared 
to find the approximate total error effect. For investigating the 
effect of measurement errors on parameter estimation accuracy, this 
method is not limited to a single-step application. For small 
errors which linearly affect the output measurements, the results 
of this simulated data analysis match those of the ensemble analysis. 


The simulated data analysis method is related to the ensemble 

analysis in that the basic equation utilized is Eq. (2.50). The 
2 

matrix 3 J is unchanged for any one step. What is different 


3P 

is that rather than computing the sensitivities 


a e 


or 

a j 


a e 


the analysis computes the residual (y. - £.) in — . 

1 1 3 p 


For output measurement biases, the residual is simply 


(y. - y i > = y 1± - y ± - B 

where y is the i^ sample of y^ defined by Equation (2.30). 
For transformation errors due to the T matrix. 


(y i " y i } = T(1L \ + Du i } “ i + Du i> = <T " I) (Hx. + Du.) 


(2.74) 


(2.75) 
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If initial conditions are unestimated, both x and x in 
Eqs . (2.1) and (2.3) are integrated, and 



(y T . ~ 9.) = «(x - x ) 

ii 1 11 


(2.76) 


For biases to the control measurements, 
(2,4) are reevaluated using 


Eqs. (2.3) and 


u . = u . + B 
mi i c 


(2.77) 


For scale factor errors, the control 


u = T u. 
mi c l 


(2.78) 


is used. 


For random errors, the total parameter covariance is 


E {6p<$p } = E 


noise 


+ I <5 V p j I) 

j=l 


(2.79) 


where r is the number of random error sources. E noise is the 
covariance of the parameter estimate errors due to the output 
noise w^. This again comes from Eq. (2.12). 
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Mean errors include some elements in the T matrix and 
all elements in and F c (Eqs. (2.32) and (2.35). If 
there are no lag errors, the evaluation is just like Eq. (2.75). 

If there are control measurement lag errors, Eq. (2.35) needs 
to be evaluated with u = u. If there are output lags, Eq. (2.32) 
needs evaluation. This requires reintegration of Eq. (2.3) 
and Eq. (2.35) in the form 

y = -F y + F ( [Hx + Du]) 

L m L m 


Equations (2.51) - (2.57) need to be reevaluated in case of 
control lags. The total effect due to lags on Eq. (2.50) is 


y i - h - - (Hi i + D “mi> 


For mean errors, the total effect on the parameter values 
is found from 


m 

E jfipf = J ( 6 Pj) 

j = 1 

where m is the number of mean error sources, and 6p are 

. j 

the mean errors due to each source computed from Eq. (2.50). 


(2.80) 


(2.81) 


(2.82) 



2.5.2 Nonlinear Analysis 

For large or non-linear measurement errors and the presence 
of process noise, a Monte Carlo simulated data analysis technique 
should be used. In this method, several different data sequences 
are simulated and used sequentially in the identification 
process. 

The random errors contained in B, part of T, B^T^, and 

x are generated at the beginning of each simulation using 

the errors' standard deviations and a random number generator. 

These errors are held constant during each single Monte Carlo 

run, but are changed from run to run. The random noise and 

w are regenerated at each sample point during each run. Each 
ci 

of the mean errors in T plus elements of and F £ are set 

equal to the constant mean values and are not changed during 
any of the runs. 

For output measurement errors in T, B and only, the 
residual (y ± - ^) in Eq. (2.50) is computed by 

' y - £ = TCHx^ + Du i ) + B + - (HSL + Du^) 


For random initial conditions, Eq. (2.1) must be integrated 
each time and Eq. (2.83) gets changed to 


- y ± = tchx. + 


Du t ) 


+ B 


+ w. - (Hx. + 
i i 


D Ui ) 


(2.83) 


(2.84) 
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For non-lag control measurement errors, Eqs. (2.34) and (2.36) 
get combined so that at each sample point 


u 


mi 


T u 
c 


+ B + 


w 


Cl 


(2.85) 


where w is randomly generated each time point. Because of 
this change, Eqs. (2.51) - (2.53) require integration each 
pass through, and Eqs. (2.54) - (2.57) require re-evaluation 
each pass through. With these changes, Eq. (2.85; becomes 


y. - y. = T (Hx. + Du.) + B + w 

11 i i i 



( 2 . 86 ) 


The error Ap in the parameter vector obtained from each run 
J 

is saved. For m Monte Carlo runs, the mean error in p is 


m 

Ap A E {Ap} = - Y Ap. 

The standard deviation about this mean is 


(2.87) 


E_ {6p6p^} 


m 2 (A P, ~ A P) ( A P, ~ Ap) T 
j=l J 


( 2 . 88 ) 
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2.6. Implementation of the Analysis Techniques 

The ensemble analysis and simulated data analyses techniques 
described in Sections 2.4 and 2.5 were coded into a digital 
computer program. This program enables the assessment of 
uncertainty (due to instrumentation errors) in the accuracy 
of the aircraft parameters identified from flight test data. 

A summary of the equations contained within this computer program 
is presented in the Appendix of this report. The longitudinal 
equations including the short period mode and the lateral equations 
are both contained in the program. 

This program has been exercised using stability and control 
derivatives from the DC-8, the F-4C, the Cessna 172, and the 
HL-10 lifting body. In all cases, an input sequence is first 
found such that the recorded output has an appropriate amount 
of information to allow the identification process to take place. 

In the next section of this report, the results of exercising 
the ensemble analysis option of the program using the F-4C as 
an example aircraft are presented. In addition to output 
measurement noise, the measurement errors which are studied 
include output biases and elements in the T matrix. 
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Ill 


STUDY OF THE EFFECT OF MEASUREMENT ERRORS ON PARAMETER IDENTIFICATION ACCURACY 

In this section, the results of a study using the previously described 
error analysis program that determines the effects of measurement errors on 
parameter identification accuracy are presented. Both longitudinal and 
lateral motion of an F-4C aircraft with typical control surface deflections 
are studied. A range of measurement and recording errors, representing the 
current state of flight instrumentation is investigated. 

3.1 Current Flight Instrumentation Accuracy 

The general measurement accuracy range of instruments for sensing air- 
craft longitudinal state outputs which is typical of current flight tests is 
presented in Table 1. A similar table representing lateral instrumentation 
accuracy is presented later as Table 10. The error sources contained in 
these tables are a result of examining product literature and the specifications 
used by flight testing organizations. Although cases were found where the 
standard deviations of instrument errors exceeded the minimum and maximum table 
values, these values are judged to be reasonable ones for this study. 

The "max" values for noise, bias, and scale factor error shown in Table 1 
are basically 0.5% of the highest dynamic range typically used in flight testing. 
Most instruments are considerably better than this; however, analog data acquisi- 
tion systems have an average accuracy of about 0.5%, so it was selected as the 
worst case. Noise and bias are related because data trace values taken before 
a maneuver is executed are used as the null points. The ability to determine 
these zero values is a direct function of both the noise and bias present. 

For a digital data acquisition system, the 0.5% error is too large. Common 
error values in the measurements due to a 10-bit data acquisition system will be 
about 0.05% and correspondingly lower for more bits. With this recording accuracy, 
the instrument errors start replacing the data acquisition errors as the important 
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Table 1 


LONGITUDINAL INSTRUMENT ERRORS 
STANDARD DEVIATIONS 


Instrument 

Units 

Full Scale 
Deflection* 

Random Noise 

Random Biases 

Random 

Scale 

Factors 

Minimum 

Base 

Maximum 


Base 

Maximum 

Minimum 

Base 

Maximum 

Pitch attitude gyro 

deg 

30-90 

.015 

.15 

.45 

.015 

.15 

.45 

.05% 

.5% 

.5% 

Pitch rate gyro 

deg/sec 

30-60 

.015 

.10 

.30 

.015 

.10 

.30 

.05% 

.5% 

.5% 

Angle-of-attack vane 

deg 

20 


.10 

.10 

.010 

.10 

.10 

.05% 

.5% 

.5% 

Pitot tube 

ft/sec 

1000 

.50 

1.00 

2.50 

.50 

1.00 

2.50 

.05% 

.5% 

.5% 

Forward accelerometer 

g's 

1 

.0005 

.005 

.005 

.0005 

.005 

.005 

.05% 

.5% 

.5% 

Vertical accelerometer 

g's 

2.5 

.0025 

.005 

.025 

.0025 

.005 

.025 

.05% 

.5% 

.5% 

Pitch accelerometer 

deg/ sec^ 

30-60 

.015 

.10 

.30 

.015 

.10 

.30 

.05% 

.5% 

.5% 


Other Types of Instrument Errors 

Minimum 

Base 

Maximum 

Instrument Lags 

Bandwidth 

(sec -1 ) 

Standard 

Deviation 

(%) 

Linear accelerometer misalignment (deg) 

.15 

.60 

.90 

Gyros 

150 

5 

Pitch angular acceleration sensitivity 
to linear acceleration (deg/sec^.g) 

.10 


.60 

Linear accelerometers 

600 

5 

Pitch gyro mass unbalance (deg/sec*g) 

.025 


.60 

Angular accelerometers 

180 


Center-of-gravity uncertainty (ft) 

.25 

.50 

1.0 

Tape recorder 

6 

5 


When a range of deflection values is given, the lower number is associated with the minimum random errors. The larger 
number is associated with the maximum random errors. 





































error sources. 


The "min" values presented in Table 1 are based on the 
lower values of References 11-13. Manufacturers' guarantees were interpreted 
as 2 a values, although in no case was the data given with any statistical 
information. 

Neither the instrument manufacturers nor the flight test agencies compile 
statistical data of instrument accuracies in the form required for this error 
analysis. The validity of the range of accuracies available, such as those 
presented in Table 1 must be questioned without supporting laboratory test 
data. 

As a means of having a reference set of instruments with which to conduct 
the study, a "baseline" set of instrument accuracies was chosen within the 
range of Table 1. This set of accuracies is listed in the "base" columns, and is 
assumed to represent values of a typical flight test program. 

The effects of instrument lags, control surface deflection measurement 

errors, angular accelerometer sensitivity to linear accelerations, and gyro 

mass unbalance were not studied in this preliminary investigation. Other 

unknown meaurement errors might exist because of voltage supply fluctuations, 

temperature effects, aircraft body bending, and nonlinear errors. Additional 

errors exist in comparing parameters obtained from wind tunnel and flight 

tests because of the uncertainty in the aircraft inertia terms which can be 

(14) 

in error up to about 5% of actual values. 

3.2 Model of the F-4C Aircraft 

To conduct the study, the F-4C aircraft was chosen with a level flight 
path and an air speed of 827 ft. /sec. The linearized perturbation equations 
about this condition, as represented by Eqs. (2.16) and (2.25) for longitudinal 
and lateral motion, were utilized. The corresponding stability and control 
derivatives of the F-4 are presented in Table 2.^^ 
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Table 2 


REFERENCE STABILITY AND CONTROL DERIVATIVE PARAMETERS FOR 

THE F - 4C AIRCRAFT 


Longitudinal 

Motion 

Lateral 
Mo t ion 

Reference 
Flight Path 

Parameter 

Nominal 

Value 

Parameter 

Nominal 

Value 

Parameter 

Nominal 

Value 

Mq 

-.719 see -1 

Y S 

-.157 sec -1 

V 

827. ft/sec 

Mw 

-.591 /ft. sec 

Y6a 

-.00338 sec ^ 

a 

0 

O 

2.6 

Mu 

O 

-.C295 /ft. sec 

Y<$r 

.0246 sec ^ 

0 

O 

2.6 

MSe 

-16.2 sec ^ 

L6 

-15.98 sec ~ 2 



Zw 

-.762 sec 1 

Lp 

-1.608 sec ^ 



Zu 

-.0617 sec 1 

Lr 

. 384 sec ^ 



Z<5e 

-1. 24ft /deg. sec 

L6a 

10.92 sec ^ 



Xw 

.0273 sec ^ 

L. 

5 r 

2.54 sec ^ 



Xu 

.00701 sec" 1 

N6 

6.563 sec ^ 





Np 

-.0997 sec ^ 





Nr 

-.343 sec ^ 





N_ 

. 707 sec ^ 





6a 






N 

6 r 

-3.902 sec^ 

- 
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The longitudinal equations of motion have the characteristic equation 

s 4 + 1.488s 3 + 9.091s 2 + .05997s + .03284 = 0 (3.1) 

Factoring this equation results in a short-period frequency of 3.01 rad/sec 

with a damping ratio of 0.246. The phugoid frequency is 0.0190 rad/sec with 
a damping ratio of 0.158. 

The lateral equations of motion have the characteristic equation 

s 4 + 2.108s 3 + 7.458s 2 + 12.86s + .1153 = 0 (3.2) 

This produces a Dutch roll frequency of 2.63 rad/sec with a damping term 
of 0.0519. The roll subsidence has a time constant of 0.548 sec. The spiral 
convergence time constant is stable with a value of 111. sec. 

3.3 Effect of Longitudinal Measurement Errors 

The reference maneuver used to identify the F-4's longitudinal stability 

and control derivative is shown in Fig. 4. The measurement data sequence 

consisted of 300 points taken every 0.05 sec for a 15 sec. time span. The 

+ ° 

elevator deflection consisted of a doublet of -2.50 followed by step inputs 

O O 

of -0.5 and 0.5 . Figure 4 shows the resulting trajectories for pitch angle, 
pitch rate, angle-of-attack, and longitudinal speed perturbations about the 
reference flight path. This sequence was selected because it provides adequate 
information for the identification process. 

3.3.1 Basic Instrumentation Error Effects 

In studying the effect of instrument errors, two different identification 
cases were used. In the first, it was assumed that only the stability and 
control derivatives were identified, so that all bias errors affected the 
total estimation uncertainty. In the second, it was assumed that -state initial 
conditions and instrument biases were estimated so that their contributions 
were essentially eliminated. In both cases, initial conditions were not used 
as error sources . 
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FIGURE 


For the full longitudinal mode, the system equations are so structured 
that the longitudinal speed perturbation's initial value appears to have almost 
the same effect on the output equations as does a bias in the pitot tube or air 
speed indicator. In other words, both of these parameters are not simultaneously 
observable from a data sequence over a short time span. Thus, for the second 
identification problem described above, the pitot tube bias is not identified, 
and it enters an error source. Most of this error is identified as a forward 
speed initial condition. 

Table 3 presents the results of using the ensemble analysis program to 
compute the standard deviations of the longitudinal parameters as identified 
from the trajectory in Fig. 4. These results are for the baseline, minimum, 
and maximum values of the random error sources listed in Table 1. The resulting 
parameter deviations are those due to white noise only and those due to the 
sum of white noise and the rest of the instrumentation errors assumed present. 

The important quantities that can be obtained from this table are the ratio of 
parameter deviation to the nominal parameter value, the increase in the deviation 
size due to error sources which aren't noise, and the effect that estimating 
biases has on the total deviations. 

From Table 3, the following conclusions can be made: 

1. Addition of non-noise error sources has a substantial effect on the 

standard deviation of the parameter estimate accuracy. As seen in 

Table 3, the errors in accuracy of parameters Mu, Zu, Xw, and Xu 
are increased by over an order of magnitude by the non-modeled 
instrument errors. This is true for the minimum, baseline > and 
maximum error values. For example, the deviation of Mu goes from 

8.7% to 119.6% (.00257 to .0353) of the parameter value for the 

baseline error set without biases being estimated. This growth 
in the standard deviations is illustrated more distinctly in the bar 
graph in Fig, 5. It must be poxnted out that the largest errors 
are in the parameters associated with the phugoid mode. This is 
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Table 3 


STANDARD DEVIATIONS OF PARAMETER ESTIMATES DUE TO INSTRUMENT ERRORS 
LONGITUDINAL EQUATIONS OF MOTION 


Parameter 

Nominal** 

Value 

Baselii 

le Errors 



Minimum Errors 

| Maximum Errors 

Blase 

Est: 

!S Not 
mated 

Bi. 

Estj 

ises* 

.mated 

Biases Not 
Estimated 

Biases* 

Estimated 

Biases Not 
Estimated 

Biases* 

Estimated 

Noise 

Only 

Total 

Errors 

Noise 

Only 

Total 

Errors 

Noise 

Only 

Total 

Errors 

Noise 

Only 

Total 

Errors 

Noise 

Only 

Total 

Errors 

Noise 

Only 

Total 

Errors 

Mq 

-.719 a -1 

. 182-2 

.704-2 

.189-2 

.683-2 

.420-3 

.271-2 

.431-3 

.265-2 

.623-2 

.162-1 

.649-2 

.128-1 

Mw 

-.591d/f-s 

.359-3 

.130-2 

.455-3 

.495-3 

.677-4 

.125-3 

.789-4 

.105-3 

.110-2 

.187-2 

.132-2 

.161-2 

Mu 

-.0295d/f *s 

.257-2 

.353-1 

.442-2 

.588-2 

.294-3 

.409-2 

.522-3 

.636-3 

.298-2 

.413-1 

.534-2 

.627-2 

M6e 

-16.2 s“ 2 

. 106-1 

.788-1 

.117-1 

.778-1 

. 166-2 

.832-2 

.183-2 

.783-2 

.306-1 

.118 

.346-1 

.810-1 

Zw 

-.762 s -i 

.914-3 

.790-2 

. 103-2 

.717-2 

.322-3 

.285-2 

.339-3 

.265-2 

.372-2 

.184-1 

.407-2 

.120-1. 

Zu 

-.0617 s _1 

.345-2 

.484-1 

.617-2 

.753-2 

.459-3 

.633-2 

.800-3 

.109-2 

.478-2 

.669-1 

.852-2 

.971-2 

Z6e 

-1.24f/d*s 

.167-1 

.150 

.183-1 

.144 

.775-2 

.734-1 

.795-2 

.715-1 

.800-1 

.365 

.843-1 

.297 

Xw 

.0273 s 1 

.494-3 

.675-2 

.769-3 

.852-2 

.605-4 

.222-2 

.777-4 

.272-2 

.588-3 

.121-1 

.779-3 

.140-1 

Xu 

.00701 s x 

.478-3 

.596-2 

.116-2 

.376-2 

.596-4 

.885-3 

.121-3 

.180-2 

.574-3 

.818-2 

. 119-2 

.727-2 


* Biases of six instruments estimated. Pitot tube bias is not directly estimated. 

** Dimensions are: s-sec; d-deg; f-ft. 

# Deviations are in the form .182-2 which means .182 x 10 -2 
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PARAMETERS 


lVIATIONS OF LONGITUDINAL MOTION 
)R BASELINE INSTRUMENT ERRORS. 






because the phugoid period is 330 sec, so the 15 sec data span 
doesn't have as much information content to obtain better accuracies 
for the parameters which govern the phugoid motion. 

2. The estimation of biases increases the parameter deviation due to 
noise only, but generally reduces the total deviation of each para- 
meter. For the baseline case, the deviations of Mw, Mu, and Zu are 
reduced more than 60%. However, by estimating biases, the deviation 
of Xw increases 26%, in the baseline case. This is because the sensiti 
vity of Xw to the misalignment in the longitudinal accelerometer 
increases when biases are estimated, and for Xw, this misalignment 

is the dominant error source. This points out that it might be better 
to structure the identification scheme so that other errors, such 
as the accelerometer misalignments, are also estimated, in addition 
to the biases. 

3. The general effects of error sources other than noise, and the effects 
of estimating instrument biases are the same for the minimum, baseline 
and maximum error values. Thus, the trends exhibited by the baseline 
error magnitude study can be used as general results. 


From Eq.(2-75) t the effect of any small random instrument error e on any 
parameter's variance can be written as 


E{6p2, total = E (5p2) noise + 


E{Sp2> other + 


(IP) « <* 2 > 


(3.3 ) 


errors 

Thus, to provide high quality parameter estimates, it is necessary to keep 
small or to keep the sensitivity of a parameter's deviation 


8e 


to an error source small. 


The sensitivities of the longitudinal parameter estimates to random and 
mean instrument errors for the baseline case when biases aren't estimated are 
presented in Table 4. Corresponding sensitivities of these parameters for the 
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Table * 


LONGITUDINAL PARAMETER SENSITIVITY TO INSTRUMENT ERRORS 
WHEN NO BIASES ARE ESTIMATED 



Random Biases 


.249-02 

.321-02 

.149-01 

.938-01 

.108+00 

.413-03 

.416-02 

.653-01 

.138-00 


.295-02 

.171-02 

,336-02 

.334-01 

.359-01 

.229-03 

.131-02 

.529-02 

.166-01 


.622-02 

-.103-01 

-. 121-02 

-.298-00 

-.395-00 

.828-02 

.368-02 

-.911-01 

-.189-00 


.258-05 

-.305-04 

.339-03 

.899-03 

.103-03 

-.250-02 

.375-02 

.410-02 

-.564-02 


.559-02 

-.294-01 

-.151-00 

-.966-00 

-.129+01 

-.967-00 

.125+00 

-.775-00 

.131+01 


.674-01 

-.852-01 

.299-00 

-. 202+01 

-.415+01 

.183-00 

-.687-01 

.758-00 

-.628+01 


Random Scale Factors 



.878-02 

.121-01 

.634-01 

.382-00 

.445-00 

.229-02 

.172-01 

.313-00 

.619-00 


.190-01 

.314-02 

.363-01 

.481-02 

.714-02 

.315-03 

.527-03 

.586-00 

.104+00 


.196-01 

.812-02 

.172-01 

.327-01 

.434-01 

.165-02 

.192-02 

.232-00 

.795-01 


.201-02 

-.639-03 

-.930-02 

-.344-01 

-. 210-01 

.632-01 

-.588-01 

-.633-01 

.122+00 


.441-03 

.223-03 

.103-02 

.110-01 

.153-01 

.778-03 

.166-01 

-.109-01 

.109-01 


.610-00 

-.167-01 

-.722-00 

.509-00 

.546-00 

-.539-01 

.376-01 

.285-00 

-.224+01 


Random Misalignments and Center-of-Gravity 
Position Errors 


Parameter 

Y nx 

Y nz 

E cgx 

e cgz 

Mq 

-.309-03 

-.372-03 

-.103-01 

-.121-03 

Mw 

-.302-04 

.605-05 

.207-05 

.700-04 

Zw 

-.124-03 

.432-03 

.103-01 

.154-03 

Mu 

-.325-02 

-.462-03 

.809-06 

.160-02 

Zu 

-.465-02 

-.678-03 

.458-03 

.223-02 

Xu 

.327-02 

.463-04 

.505-05 

-.505-03 

Xw 

.882-02 

-.264-04 

.252-04 

-.287-02 

M<$e 

-.389-02 

-.203-03 

.229-03 

.899-03 

Z<5e 

.209-02 

.756-03 

.284-00 

-.242-02 


Mean Errors — Accelerometer Positions 
and Angle-of-Attack Vane Position 


Parameter 

e ax 

e az 

fc vx 

Mq 

-.103-01 

-.121-03 

.154-04 

Mw 

.624-05 

.700-04 

-.417-05 

Zw 

.103-01 

.154-03 

-.294-05 

Mu 

.144-03 

.160-02 

-.143-03 

Zu 

.648-03 

.223-02 

-.189-03 

Xu 

.687-06 

-.505-03 

.436-05 

Xw 

.224-04 

-.287-02 

.276-05 

M Se 

.675-04 

.899-03 

; .162-03 

Z <5e 

.284-00 

-.242-02 

.115-03 
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baseline case in which biases are estimated are presented in Table 5. By 
using Eq. (3.3) and the range of non-noise errors of Table 1, one can determine 
which error sources have the major effect on the accuracy of each parameter 
estimated. This was done, and the results are delineated in Table 6. 

It is noted from Tables 4, 5, and 6 that adding the capability of estimating 
initial conditions and biases tends to restructure the values of the elements in 
the sensitivity matrix. As shown in Table 6, directly estimating all the instru- 
ment biases except the pitot tube bias b^ causes b to emerge as a 
dominant error in the estimation of Mw, Mu, and Zu. But, the same estimation 
scheme removes b^ as a major error source of Xw and Xu. 

Sensitivity tables such as Tables 4 and 5 are useful in specifying the 
accuracy required of the instruments or other aircraft parameters which affect 
the accuracy of the estimated stability and control derivatives. As an example 
of this application. Fig. 6 illustrates the deviation of the parameter Z6e due 
to the uncertainty in the longitudinal position of the aircraft center-of-gravity . 
For the ZSe uncertainty to be less than 10% of the nominal value, the position 
of the center-of-gravity must be known to within 0.4 ft. 

3-3.2 Effect of Changed Input and Data Span 

There are several other effects which must be considered in drawing general 
conclusions of the importance of instrumentation errors on flight identification. 
Some of these include : 

The type and configuration of the aircraft being studied 
The control input sequence used to excite the aircraft 
The sampling rate and time span of the data collected 
The types of instruments available 

Doubling the amplitude of elevator deflection cuts the effect of bias errors 
in half. All other error effects are unchanged by the increased amplitude. This 


1 . 

2 . 

3. 

4. 
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Table 5 


LONGITUDINAL PARAMETER SENSITIVITY TO INSTRUMENT ERRORS WHEN 
INITIAL CONDITIONS AND BIASES EXCEPT b u ARE ESTIMATED 


Random Scale Factors 

Parameter 

e e 


e a 

e u 

e nx 

e nz 

e q 

Mq 

.486-02 

-.107-02 

-.254-01 

-.248-02 

-.113-02 

.585-00 

-.560-00 

Mw 

-.469-03 

.441-02 

.229-01 

.195-02 

.165-03 

-.290-02 

-.261-01 

Zw 

-.198-01 

.456-01 

.229-01 

.220-02 

.734-03 

-.712-00 

.661-00 

Mu 

-.134-00 

-.454-01 

.385-00 

.480-01 

.211-02 

.260-00 

-.515-00 

Zu 

-.142-00 

-.811-01 

.529-00 

.899-01 

.143-02 

.171-00 

-.569-00 

Xu 

.124-02 

.326-02 

-.184-01 

.933-02 

-.871-02 

-.101-01 

.234-01 

Xw 

.563-03 

-.187-02 

-.691-03 

-.171-02 

.266-01 

.223-02 

-.252-01 

M<5e 

-.204-00 

-.510-00 

-.179-00 

.228-02 

.135-02 

.486-01 

-.154+02 

Z<5e 



.186-00 

.741-01 

.787-01 

.259-01 

-.440-02 

-.239+01 

.783-00 


Random Bias, Misalignment 

s, and Center-of-Gravity 



Position 

Errors 



Parameter 

b u 

^nx 

*r- 

'nz 

„ 

Cgx 

e cgz 

Mq 

-.137-03 


-.352-03 

-.103-01 

.172-03 

Mw 

.557-04 

.891-04 

-.296-05 

.109-04 

.646-04 

Zw 

.249-03 

.453-03 

.423-03 

.103-01 

.503-04 

Mu 

.137-02 

.120-02 

-.281-03 

-.212-04 

-.409-03 

Zu 

.129-02 


-.384-03 

.438-03 

-.387-03 

Xu 


-.715-03 

.125-04 

.150-04 

.710-02 

Xw 

-.132-04 


-.516-07 

.190-04 

-.739-02 

M6e 

.784-03 

.779-03 

-.749-04 

.613-04 

-.164-02 

Z6e 

-.196-02 

-.308-02 

.893-03 

.284-00 

-.242-02 


Mean Errors — Accelerometer Positions 
and Angle-of-Attack Vane Position 

Parameter 

£ ax 

e az 

e vx 

Mq 

-.103-01 

.172-03 

.172-04 

Mw 

.354-04 

.646-04 

-.245-04 

Zw 

.103-01 

.503-04 

-.146-04 

Mu 

.502-03 

-.409-03 

-.524-03 

Zu 

.116-02 

-.387-03 

-.722-03 

Xu 

-.971-05 

.710-02 

.247-04 

Xw 

.193-04 

-.739-02 

-.293-06 

M6e 

-.130-03 

-.164-02 

.191-03 

Z6e 

.284-00 

-.242-02 

.506-04 
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Table 6 


DOMINANT RANDOM ERROR SOURCES FOR LONGITUDINAL PARAMETER IDENTIFICATION 


Parameter 

Biases . Not 
Estimated 

All Biases 
but bu 
Estimated 

Mq 

e , e . , e 
nz q cgx 

e nz ’ e q ’ £ cgx 

Mw 

b 

a 

b u ’ e a* e q 

Mu 

b 

a 

noise, b u 

M6e 

e q 

e q 

Zw 

b 0’ £ cgx 

e nz* G cgx 

Zu 

b 6‘ b a 

noise, b u 

Z6e 

e xcg 

£ xcg 

Xw 

b u> \x 

Y 

'nx 

Xu 

b u 

e 

cgz 
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STANDARD DEVIATION OF Z6e/ INOMINAL VALUE 






can be directly obtained from Eq. (2.59). 


Another effect of input changes was found by doubling the pulse length 
of the elevator deflections shown in Fig. 4 but halving the pulse amplitude. 

The data sequence was doubled to 30 sec with samples taken 0.10 sec. The 
resulting input had the same number of sample points and same area under the 
input deflection curve. The results of the input time history as compared to 
the baseline case are illustrated in Table 7. When biases are not estimated, 
doubling the time span increased the noise-only and total errors for 4 para- 
meters; it decreases the noise-only and total errors for three parameters; 
and it increases the noise effect but reduces the total errors for two of the 
parameters. Notice that phugoid mode parameter accuracy is generally improved 
by increasing the time span, and the accuracy of the short period mode para- 
meters is generally reduced. 

As shown in Table 7>, when biases are estimated, doubling the time span 
and halving the input amplitude decreases the accuracy of all parameter estimates. 

For this case, the combined effect of the reduced signal-to-noise ratio and the 
sensitivity changes due to the addition of more parameters being estimated resulted 
in no improved accuracies. Table 7 points out that any definite conclusions on 
instrument accuracy effects are dependent upon the maneuvers flown and the para- 
meters being estimated. 

3.3.3 Effect of Fewer Instruments 

The above results were obtained assuming that the aircraft has seven instru- 
ments for obtaining longitudinal information. It was assumed that these instruments, 
discussed in Section 2.3, had errors as modeled in Eqs . (2. 60)-(2. 61) . Parameter 
identification can be conducted with fewer instruments, and so it is desirable 
to know what reducing the number of instruments has on the overall parameter 
accuracy. It is known that reduced instrumentation increases the parameter 
uncertainty due to white noise only. But it is conceivable that removing an 
instrument also removes a major unmodeled and unestimated error source. 


To test this idea, a set of runs was made in which different instruments 



Table 7 


EFFECT OF DOUBLING THE TIME SPAN OF THE MEASUREMENT SEQUENCE WITH CONSTANT INPUT ENERGY 



Reference Trajectory 

Doubled Time Span 


Biases Not 

Biases 


Biases I 

Mot 

Biases 

Nominal Value 

Estimated 

Estimated 


Estimate! 

i 

Estimated 

L i 


Noise 

Total 

Noise 

Total 

Noise 

Total 

Noise 

Total 


Only 

Errors 

Only 

Errors 

Only 

Errors 

Only 

Errors 

-.719 s -1 

.182-2 

.704-2 

.189-2 

.683-2 

.282-2 

.787-2 

.333-2 

.724-2 

-.591d/f-s 

.359-3 

. 130-2 

.455-3 

.495-3 ! 

.432-3 

.784-3 

.123-2 

.128-2 

-.0295d/f *s 

.257-2 

. 353-1 

.442-2 

.588-2 

.158-2 

.214-1 

.255-1 

.270-1 

1 

-16.2s" 2 

. .106-1 

.788-1 

.117-1 

.778-1 

.221-1 

.837-1 

.307-1 

.798-1 

-.762 s -1 

. 914-3 

.790-2 

.103-2 

.717-2 

.142-2 

.784-2 

.211-2 

.718-2 

-.0617s -1 

.345-2 

.484-1 

.617-2 

.753-2 

.213-2 

.295-1 

.336-1 

.353-1 

-1.24f/d*s 2 

.167-1 

.150 

.183-1 

.144 

.342-1 

.191 

.482-1 

.153 

.0273s" 1 

.494-3 

.675-2 

.769-3 

.852-2 

.698-3 

.739-2 

.128-2- 

.938-2 

.00701 s -1 

.478-3 

.596-2 

.116-2 

.376-2 

.238-3 

.249-2 

.581-2 

.597-2 














were individually and collectively removed from use. The resulting deviation 
of the parameters is shown in Table 8 for cases with bias estimated and not 
estimated. The comments which can be made from Table 8 are: 

1. Removal of the angle-of-attack vane approximately doubles the deviation 
of Mu and Zu for both cases where biases are and are not estimated. 

This is substantial m terms of the nominal values of these parameters. 

The percentage change in the deviation of Mw is large, but the values 
of the deviations are small with respect to the nominal value (.591). 

2. Removal of the pitch angular accelerometer alone does not substantially 
affect any of the parameter deviations. However, flight test personnel 
have commented that they only correct for center-of -gravity position 
errors in other inertial instruments when measurements of the angular 
accelerations are directly available. So, in a practical sense, no 
pitch accelerometer also means additional errors elsewhere. 

3. Removal of the pitot tube alone when biases are not estimated has a 
large effect on the total deviations of Xw and Xu. The removal increases 
the deviation percentages from 24.7% and 85.1% to 32.8% and 125.8% for 

Xw and Xu, respectively. Removal of the pitot tube when biases are esti- 
mated has the largest effect on Mu, Zu, and Xu. However, the deviations 
of these parameters with biases estimated are all smaller than when the 
pitot tube data are used without estimating the biases. 

4. Removal of the angle-of-attack vane, the pitot tube, and the pitch angular 
accelerometer simultaneously makes it highly desirable to estimate biases. 
For this situation, the parameters Mu, Zu, and Xu again were most highly 
affected by the removal of the three instruments. The deviations of 

all parameters but Zw and Z6e were increased by the instrument removal. 

The data shown in Tables 3, 4, 5, and 8 are quite useful in specifying what 
the best set of instruments are and what the corresponding instrument accuracies 
must be to obtain parameter accuracies within some acceptable level. Figure 7 
illustrates the effect of instrument bias variations on the accuracy of Mu for 
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Table 8 


EFFECT OF INSTRUMENT REMOVAL ON FARAMETER DEVIATIONS FOR BASELINE SET 
OF ERRORS AND LONGITUDINAL EQUATIONS OF MOTION 


Parameter 

** 

Nominal 

Value 

All 

Instruments 

No Angle-of- 
Attack Vane 
(a) 

No Pitch Angle 
Accelerometer 
(q) 

No Pitot Tube 

m 

No a 
or 

» q, 

u 

Biases 

No 

Biases 

Biases 

No 

Biases 

Biases 

bu 

Biases 

bu 

Biases 

bu 

Mq 

-. 719s _1 

.704-2 

.683-2 

.766-2 

.679-2 

.120-1 

.811-2 

.705-2 

.676-2 

.108-1 

.825-2 

Mw 

- . 591d/ f • s 

.130-2 

.495-3 

.291-2 

.690-3 

.124-2 

.669-3 

.130-2 

.869-3 

.290-2 

.171-2 

Mu 

- . 0295d/ f • s 

O 

.353-1 

.588-2 

.679-1 

.138-1 

•361-1 

.589-2 

.354-1 

.182-1 

.716-1 

.400-1 

M<5e 

-16.2s 

1 

.788-1 

.778-1 

.792-1 

.785-1 

.181 

.837-1 

.788-1 

.784-1 

.151 

.953-1 

Zw 

-.762s 

.790-2 

.717-2 

.815-2 

.717-2 

.926-2 

.738-2 

.792-2 

.712-2 

.780-2 

.755-2 

Zu 

-.0617s 

.484-1 

.753-2 

.921-1 

.181-1 

.490-1 

.763-2 

.484-1 

.318-1 

.960-1 

.527-1 

Z6e 

-1.24f/d-s 

.150 

.144 

.150 

.143 

.148 

.144 

.150 

.144 

.148 

.144 

Xw 

.0273s 

.675-2 

.852-2 

.681-2 

.852-2 

.674-2 

.852-2 

.896-2 

.878-2 

.897-2 

.878-2 

Xu 

.00701s 

.596-2 

.376-2 

.608-2 

.384-2 

.596-2 

.376-2 

.882-2 

.505-2 

.937-2 

.658-2 


Biases" means that all biases are error sources; "bu" means all biases except b are directly estimated. 
** u 

s = second; d = degrees; f = feet. 


Dimensions are: 




STANDARD DEVIATION OF Mu / NOMINAL VALUE 



0 min .1 .2 .3 .4 max .5 

PITCH ATTITUDE GYRO BIAS - DEG 

III. I I II. 

0 min .02 .04 .06 .08 .10 

ANGLE-OF-ATTACK VANE BIAS - DEG max 

lit a i i. I 

0 min .005 .010 .015 .020 .025 

max 

VERTICAL ACCELEROMETER BIAS - g's 


FIGURE 7 EFFECT OF VARIOUS INSTRUMENT BIASES ON THE ESTIMATION ACCURACY 
OF THE PARAMETER Mu. ALSO SHOWN IS THE EFFECT OF NOT USING 
THE ANGLE-OF-ATTACK VANE DATA. 


c 
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the baseline set of errors. The standard deviation of the vertical accelero- 
meter, pitch attitude gyro, and angle-of-attack vane biases are indicated by 
the values where their curves cross the reference error line. Reduction of 
the angle-of-attack vane bias can reduce the ratio of the standard deviation 
of Mu to its nominal value from about 1.2 to 0.65. On the other hand, an 
increase in either the vertical accelerometer or pitch attitude gyro bias 
from the baseline values can cause large increases in the error in Mu. With 
baseline error values , the removal of the angle-of-attack vane causes the 
standard deviation of Mu to almost double (ratio goes from 1.2 to 2.3). If 
biases are estimated, the standard deviation of Mu is reduced about 85%. 

Removing the use of the a-vane when biases are estimated still is considerably 
better than the reference baseline case. 

For some stability and control derivatives , it is possible to improve their 
accuracy by not using an instrument’s data. Figure 8 illustrates such a case. 
Here, the ratio of the standard deviation of Xw to its nominal value is shown 
as a function of the standard deviation of the pitot tube bias. For the base- 
line reference case, the pitot tube's deviation is one (1) ft/sec and it is 
better to use the pitot tube data. If the maximum value of the pitot tube 
bias is expected, however, (Refer to Table 1), it is better not to use this 
data in estimating Xw. The cross-over deviation of b^ beyond which the pitot 
tube shouldn't be used is about 1.8 ft/sec. 

3.3.4 Effect of Changing the Algorithm Weighting Matrix 

Referring back to Table 8 again, to the case where the three instruments 

are removed and the biases are estimated (last column), it can be seen that the 

deviations of Mq, Mw, M6e, and Zw are acceptably small. The main contributions 

to errors in Mu and Zu are due to noise. The chief error sources affecting 

Z<5e, Xw, and Xu are the center-of-gravity uncertainty and the misalignment of 
the forward accelerometer. A flight test requirement might be to improve the 
accuracies of the two gyros and two linear accelerometers so these error sources 
are acceptably reduced. This may also include better calibration of the center- 
of-gravity position. 
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STANDARD DEVIATION OF Xw / [NOMINAL VALUE 


FIGURE 8 



PITOT TUBE BIAS - b (ft / sec) 

u 


EFFECT OF PITOT TUBE BIAS ON THE ESTIMATION ACCURACY OF THE PARAMETER Xw 
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More computer runs were made to determine the effect of reducing instrument 
error values on the overall parameter accuracy when only four instruments are 
used. The results are presented in Table 9. Concentration was placed on the 
five parameters (Mu, Zu, Z<Se, Xw, and Xu) with sizeable deviations. Case (a) 
is the reference case which is the same as the last column of Table 8. In 
Table 9, the parameter deviations due to noise only, to other errors, and their 
root-sum-square total is shown. 

In Case (b) , the standard deviations of the white noise errors were halved. 
(The other error statistics were held constant.) This reduced the diagonal 
elements of the matrix R used in the cost function J (see Eq.(2.6 )) to one- 
fourth their reference value. This caused the anticipated result - the errors 
in the five parameters due to noise were cut in half, and this significantly 
affected the deviations of Mu and Zu. 

In Case (c) , the standard deviations of the four noise sources and other 
errors were set to the minimum values shown in Table 1. As was expected, the 
deviations of the parameters due to noise only were all reduced. But surprise- 
ingly , the error in Xu due to other errors increases. This is due to the fact 
that in reducing the noise values, the ratios of the elements of the weighting 
matrix R are changed. This causes the elements of the sensitivity matrix 
of parameter deviations to error sources to change. In the case of Xu, the 
sensitivity of Xu to the error in the vertical position of the center-of- 
gravity increased from 0. 0113/f fsec to 0. 0318/ft • sec . Thus, the deviation of Xu 
also grew. 

The best set of noise parameters (and cost function weighting terms) was 
found to be setting the longitudinal accelerometer noise to the baseline value 
and the other three terms to the minimum value. Results of this situation .are 
shown as Case (d) . 


The above four cases illustrate that the effect of instrumentation error 
sources on stability and control derivative estimate accuracies is also highly 
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Table 9 


EFFECT OF CHANGED NOISE VARIANCES ON TOTAL PARAMETER ESTIMATE DEVIATIONS 


M 


Parameter 

Nominal 

Value 

Reference Case 
(a)* 

Noise Parameters 
Halved 
(b)* 

Minimum Noise With 
Reduced Errors 
(c)* 

Special Case 
(d)* 

Parameter Deviations 

Parameter Deviations 

Parameter Deviations 

Parameter Deviations 

Noise 

Only 

Other 

Errors 

Total 

Noise 

Only 

Other 

Errors 

Total 

Noise 

Only 

Other 

Errors 

Total 

Noise 

Only 

Other 

Errors 

Total 

Mu 

- . 0295d/ f • s 

.369-1 

.153-1 

.401-1 

.185-1 

.153-1 

.240-1 

.528-2 

.140-1 

.149-1 

.529-2 

.136-1 

.146-1 

Zu 

-.0617s -1 

.483-1 

.212-1 

.527-1 

.242-1 

.212-1 

.321-1 

.744-2 

.235-1 

.246-1 

.811-2 

.146-1 

.167-1 

Z6e 

-1.24f/d-s 2 

.217-1 

.142 

.144 

.109-1 

.142 

.142 

.928-2 

.720-1 

.726-1 

.953-2 

.717-1 

.723-1 

Xw 

.0273s" 1 

.783-3 

.875-2 

.878-2 

.392-3 

.875-2 

.876-2 

.793-4 

.278-2 

.278-2 

.776-3 

.273-2 

.284-2 

Xu 

.00701s" 1 

.326-2 

.571-2 

.658-2 

.163-2 

.571-2 

.594-2 

.679-3 

.804-2 

.807-2 

.164-2 

.358-2 

.393-2 


(a) 

w 0 = 

.15° 

(b) 

w e 

.075 

(c) 

w e = 

.015 

(d) 

w e 

.015 


.10°/sec; w = .005 g's; w = .005 g's; y = .6°; e = .5 ft 

ax nz nx,z cgx,z 


w„ = 

. 05° / sec ; 

w 

q 


nx 

w = 

. 015° / sec 

; w 

q 

nx 


O 


w = 

.015 /sec 

; W 

q 


nx 


nz 

w 

nz 

v. = 
nz 


y = .6°; e. = .5 ft 
nx,z cgx,z 

y = .15°; e = .25 ft 

nx,z cgx,z 

Y nx,z ,1 “* » e cgx,z “ * ^ 




dependent upon the weighting matrix used by the Newton-Raphson identification 
algorithm. In fact, an important problem which should be solved is how to 
choose the weighting matrix for the identification process so that the total 
parameter deviations due to all types of measurement errors are minimized. 

3.4 Effect of Lateral Measurement Errors 

Similar program runs to those discussed in Section 3.3 were made for the 
lateral equations of motion of the F-4 aircraft. Again, seven instruments 
were modeled which are described by Eqs .(2 . 60) — (2 . 61) . . The range of instrument 
and calibration errors investigated are similar to those of the longitudinal 
instruments, and they are presented in Table 10. 

The aileron and rudder deflections used to excite the lateral motion and 
the resulting aircraft trajectory are depicted in Fig. 9. Again, a 15 sec time 
span was simulated, and 300 data points securing every 0.05 sec. were processed. 
The trajectory shown in Fig. 9 was used as the reference case for the lateral 
study. 

A study was made of the identification process applied to the reference 
trajectory using the baseline, minimum, and maximum sets of instrument errors 
listed in Table 10. The results are presented in Table 11. The conclusions 
which can be made from the lateral study with all seven instruments are: 

1. The addition of non-noise error sources to the output measurements 
causes a substantial effect on the accuracy of the estimated parameters. 
The standard deviations of Y0, Y6a, Y6r, Lp, and Np increased by over 
an order of magnitude for the baseline case. The standard deviation 

of Y6a increased from about 3% to 164% of the nominal value. These 
increases are depicted more clearly in Fig. 10. The trends are the 
same for the minimum and maximum error value cases. 

2. Unlike the longitudinal equations of motion, the addition of bias 
estimation was more detrimental than good. Of the 13 parameters 
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Table 10 


LATERAL INSTRUMENT ERRORS 
STANDARD DEVIATIONS 




Full Scale 

Random Noise 

Random Biases 

Random 

Scale 

Factors 

Instrument 

Units 

Deflection* 

Minimum 

Base 

Maximum 

Minimum 

Base 

Maximum 

Minimum 

Base 

Maximum 

Angle-of-sideslip 

vane 

deg 

10 

.005 

.05 

.05 

.005 

.05 

.05 

.05% 

.5% 

.5% 

Roll rate gyro 

deg/ sec 

30-150 

.015 

.10 

.75 

.015 

.10 

.75 

.05% 

.5% 

.5% 

Yaw rate gyro 

deg/ sec 

30-60 

.015 

.10 

.30 

.015 

.10 

.30 

.05% 

.5% 

.5% 

Roll attitude gyro 

deg 

360 

.09 

.50 

.90 

.09 

.50 

.90 

.05% 

.5% 

.5% 

Lateral accelerometer 

g's 

.2 

.0001 

.0005 

.001 

.0001 

.0005 

.001 

.05% 

.5% 

.5% 

Roll accelerometer 

deg/sec^ 

30-60 

.015 

.10 

.30 

.015 

.10 

.30 

.05% 

.5% 

.5% 

Yaw accelerometer 

deg/sec 

30-60 

.015 

.10 

.30 

.015 

.10 

.30 

.05% 

.5% 

.5% 


Other TvDes of Instrument Errors 

Minimum 

Base 

Maximum 

Instrument Lags 

Bandwidth 

(sec - -*-) 

Standard 

Deviation 

(%) 

Rate gyro misalignment (deg) 

.15 

.60 

.90 

Gyros 

150 

5 

Rate gyro mass unbalance (deg/sec’g) 

.025 


.60 

Linear accelerometers 

600 

5 

Angular accelerometer misalignment (deg) 

.15 

.60 

.90 

Angular accelerometers 

180 


Angular accelerometer sensitivity to 







linear acceleration (deg/sec^.g) 

.10 


.60 

Tape recorder 

6 

5 

Center-of-gravity uncertainty (ft) 

.25 

.5 

1.0 





*When a range of deflection values is given, the lower number is associated with the minimum random errors. The larger 
number is associated with the maximum random errors. 










Table 11 


STANDARD DEVIATIONS OF PARAMETER ESTIMATES DUE TO INSTRUMENT ERRORS 

LATERAL EQUATIONS OF MOTION 





Baseline 

Errors 



Minimum 

Errors 



Maximum 

Errors 



** 

Biases Not 
Estimated 

Biased 

Estimated 

Biases Not 
Estimated 

Biased 

Estimated 

Biases Not 
Estimated 

Biases 11 ' 

Estimated 


Nominal 

Noise 

Total 

Noise 

Total 

Noise 

Total 

Noise 

Total 

Noise 

Total 

Noise 

Total 

Parameter 

Value 

Only 

Errors 

Only 

Errors 

Only 

Errors 

Only 

Errors 

Only 

Errors 

Only 

Errors 

YB 

- . 157s _1 

.347-3 

.683-2 

.356-3 

.636-2 

.581-4 

.310-2 

.599-4 

.296-2 

.820-3 

.133-1 

.845-3 

.117-1 

Y6a 

-.00338s" 1 

.102-3 

.554-2 

.136-3 

.694-2 

.199-4 

.278-2 

.246-4 

.327-2 

.211-3 

.110-1 

.297-3 

.141-1 

Y6r 

.0246s" 1 

.117-3 

.295-2 

.124-3 

.305-2 

.230-4 

.150-2 

.242-4 

.156-2 

.239-3 

.574-2 

.254-3 

.604-2 

LB 

-15.98s -2 

.388-1 

.247 

.401-1 

.271 

.584-2 

.572-1 

.611-2 

.127 

.108 

.537 

.112 

.536 

Lp 

-1.608s -1 

.298-2 

.375-1 

.317-2 

.585-1 

.464-3 

.143-1 

.510-3 

.237-1 

.923-2 

.107 

.945-2 

.131 

Lr 

.384s -1 

.114-1 

.389-1 

.115-1 

.595-1 

.175-2 

.132-1 

.179-2 

.303-1 

.339-1 

.106 

.341-1 

.110 

LSa 

10.92s" 2 

.130-1 

.941-1 

.152-1 

.764-1 

.201-2 

.256-1 

.235-2 

.179-1 

.386-1 

.243 

.431-1 

.950-1 

i L6r 

2.54s -2 

.118-1 

.862-1 

.123-1 

.557-1 

. 180-2 

.217-1 

.188-2 

.608-1 

.345-1 

.217 

.356-1 

.144 

NB 

6.563s -2 

.341-2 

.232-1 

.370-2 

.206-1 

.553-3 

.536-2 

.598-3 

.741-2 

.872-2 

.681-1 

.964-2 

.244-1 

Np 

-.0997s" 1 

.364 -3 

.490-2 

.588-3 

.132-1 

.602-4 

.231-2 

.978-4 

.546-2 

.961-3 

.899-2 

.150-2 

.280-1 

Nr 

-.343s -1 

.143-2 

.819-2 

.166-2 

.206-1 

.234-3 

.331-2 

.271-3 

.842-2 

.371-2 

.162-1 

.432-2 

.484-1 

N6a 

.707s -2 

.144-2 

.833-2 

.234-2 

.280-1 

.232-3 

.159-2 

.383-3 

.115-1 

.379-2 

.282-1 

.598-2 

.624-1 

N6r 

-3.902s -2 

.705-2 

.497-1 

.713-2 

.434-1 

.104-2 

.113-1 

.105-2 

.956-2 

.183-1 

.752-1 

.187-1 

.647-1 


•k 

Biases of six instruments estimated. Roll angle gyro bias is not directly estimated. 

** 

Dimension is: s = second. 




RATIO OF PARAMETER UNCERTAINTY DUE TO ALL ERRORS TO UNCERTAINTY 

DUE TO NOISE 





estimated, only 4, 2, and 6 total standard deviations were lowered 
by estimating biases for the baseline, minimum, and maximum error 
sets, respectively. This is because the predominant error sources 
are not due to biases for the lateral equations. By estimating biases, 
the sensitivity of the parameters to the other error sources is generally 
increased. Thus, the resulting estimate error increases. It is con- 
cluded that for the range of instrument errors given and the trajectory 
flown it is better not to estimate biases. However, the dependence 
on the aircraft parameters and trajectory should be taken into account 
before any general statement can be made concerning whether or not 
biases should be estimated. 

Tables 12 and 13 present the aircraft lateral parameter sensitivities to 
the random and mean error sources which affect the accuracy of the system output 
measurements. As with the longitudinal sensitivities presented in Tables 4 and 
5, these sensitivities are useful for determining what the primary sources of 
error are, which affect the estimation accuracy of the lateral parameters. The 
sensitivities are also useful for specifying instrumentation and calibration 
accuracy required to provide a given level of identification accuracy. 


Table 14 presents the dominant random error sources which affect the 
lateral parameters variances. These sources were determined from use of 
Tables 10, 12, and 13. To improve the estimation accuracy of any one given 
parameter requires concentrating on lowering the magnitude of the chief error 
sources affecting that parameter. For example, for a more accurate parameter 

Np, the uncertainty in the center-of-gravity position of the aircraft should 
be reduced. 

The effect of reducing the number of instruments used to measure lateral 
motion was also studied. Table 15 contains the results of simulating the 
identification process while removing individually and collectively the angle- 
of -sideslip vane, the roll angle accelerometer, and the yaw angle accelero- 
meter. Conclusions which can be made from Table 15 are: 
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Table 12 


LATERAL PARAMETER SENSITIVITY TO INSTRUMENT ERRORS 
WHEN NO BIASES- ARE ESTIMATED 


Random Biases 


Parameter 


Y2 

L6 

NB 

Lp 

Np 

Lr 

Nr 

Y<5a 

L5a 

NS a 

Y6r 

LSr 

N5r 


.560-02 

-.160-02 

.580-00 

-.959-00 

.994-02 

-.810-01 

.112-02 

-.294-01 

.204-02 

-.835-02 

.396-03 

-.475-02 

.698-02 

-.269-01 

.166-03 

-.838-03 

.257-01 

-.255-00 

.790-02 

.923-02 

.295-03 

-.509-03 

.522-01 

-.394-00 

.137-00 

-.421-01 


-.420-03 

-.265-03 

-.231+01 

-.142-03 

-.342-04 

-.179-00 

-.143-00 

-.499+01 

-.306-01 

-.306-01 

-.163-01 

-.140-01 

-.1501-01 

-.889-02 

.201-02 

-.486-02 

-.352-02 

.275+01 

.355-02 

-.504-03 

-.189-02 

-.147-02 

-.111+01 

-.599-03 

-.215-03 

.128-01 

.153-01 

.141+02 

.163-01 

-.289-02 

-.525-02 

-.322-02 

-.136+01 

-.470-02 

.159-04 

-.241-03 

-.148-03 

-.910-00 

-.779-04 

-.445-04 

-.633-01 

-.532-01 

-.240+02 

-.295-01 

-.106-01 

.182-02 

.916-03 

.168+01 

.306-04 

.192-03 

-.183-03 

-.979-04 

-.120+01 

.201-05 

-.720-04 

-.769-01 

-.611-01 

-.438+01 

-.300-01 

-.205-01 

-.106-01 

-.622-02 

-.263+01 

-.812-02 

-.195-02 


Random Scale Factors 


e 


.720-02 

.697-00 

.214-01 

.359-02 

.213-02 

.105-01 

.398-02 

.168-03 

.272-01 

.934-02 

.262-03 

.358-01 

.171-00 


.990-03 
.240+01 
.779-01 
. 121+00 
.249-01 
.172-00 
.248-01 
.283-03 
.176+01 
.726-01 
.327-03 
.864-00 
.120+00 


.947-02 

.957-00 

.431-01 

-.886-03 

-.361-02 

-.159-01 

-.514-02 

.197-03 

.215-01 

.109-01 

.530-04 

.857-02 

-.230-00 


-.258-03 
-.144-00 
-.549-02 
.949-02 
.434-03 
-.104-01 
-.495-02 
-.980-04 
.524-02 
. 302-03 
-.766-04 
-.105+00 
-.985-02 


.157-00 

.277-01 

.362-01 

.241-02 

.658-02 

.110+00 

.152-00 

.378-02 

.638-02 

.313-02 

.246-01 

.151-01 

.736-05 


.106-01 
. 122+02 
.549-00 
.130-00 
.529-01 
.570-00 
.340-01 
.326-04 
. 908+01 
.410-00 
.159-03 
.353+01 
.285-01 


.129-00 

.131+02 

.604-00 

.671-02 

-.659-01 

-.251-00 

-.147-00 

.382-03 

.237-01 

.201-00 

-.423-03 

-.461-01 

-.334+01 


Random Misalignments an d Center-of-Gravity Position Errors 

Parameter i * ^ ’ 1 7-7 7 ^ 7 77 7 p 



.582-05 

-.293-03 

.102-03 

-.249-02 

-.967-02 

.859-02 

-.759-02 

-.396-01 

-.970-01 

-.283-00 

.109+00 

.109-01 

-.377-03 

-.221-02 

-.425-02 

-.187-01 

.319-01 

.119-02 

.401-03 

-.130-02 

.133-02 

-.685-02 

.738-01 

.393-02 

.584-05 

.224-03 

.526-05 

.179-02 

-.913-02 

-.224-03 

-.162-02 

-.214-02 

-.201-01 

.111-02 

.603-01 

.287-01 

-.198-03 

.282-03 

-.234-02 

.545-02 

-975-02 

-.812-02 

.129-05 

.108-04 

-.104-05 

.541-04 

.110-01 

.580-03 

-.238-02 

.490-02 

-.103-01 

.268-01 

-.136-00 

-.125-01 

-.233-03 

-.181-03 

-.136-03 

-.146-02 

-.154-01 

-.944-03 

-.241-05 

.104-04 

.142-05 

.116-03 

.339-02 

-.467-02 

.248-02 

-.958-03 

.652-01 

.131-01 

-.106+00 

-.529-02 

-.358-04 

.758-02 

.297-03 

.696-01 

.347-01 

.195-02 


Mean Errors — Accelerometer Positions 
and Angle-of-Sideslip Vane Position 


Parameter 



e ax 

e az 

E vx 

-.967-02 

-.859-02 

.199-05 

.109+00 

-.109-01 

.121-03 

.318-01 

-.119-02 

.443-04 

.738-01 

-.393-02 

-.826-05 

-.913-02 

.224-03 

.726-06 

.604-01 

-.287-01 

-.850-04 

-.974-02 

.812-02 

-.102-04 

.110-01 

-.580-03 

-.721-08 

-.136-00 

.125-01 

.318-04 

-.154-01 

.944-03 

.842-05 

.339-02 

.467-02 

. 260-06 

-.106+00 

.529-02 

-.272-04 

.347-01 

-.195-02 

-.689-04 
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Table 13 


LATERAL PARAMETER SENSITIVITY TO INSTRUMENT ERRORS WHEN 

INITIAL CONDITIONS AND BIASES EXCEPT b . ARE ESTIMATED 

<P 


Random Scale Factors 

Parameter 

e B 


e r 

masmm 


fP 

ef 

yb 

.671-02 

.126-02 

.959-02 


-.158-00 


.130-00 

lb 

.626-00 

-.197+01 

.958-00 


-.472-01 

-.130+02 

. 134+02 

nb 

.285-01 

-.185-01 

.460-01 


-.358-01 


.547-00 

LP 

-.674-02 

.187-00 

-.346-02 


-.271-02 


.457-01 

Np 

.766-03 

.433-01 

-.294-02 




-.865-01 

Lr 

-.119-01 

.184-02 

-.224-01 


-.106+00 


-.212-00 

Nr 

.159-02 

.147-01 

-.392-02 


.153-00 


-.184-00 

Y6a 

-.352-03 

-.505-02 

.256-03 


-.407-02 


.371-02 

Lfia 

-.194-01 

.132+01 

.337-01 

.514-01 

-.195-01 

.926+01 

.292-00 

N<Sa 

.492-02 

.168-00 

.480-02 

.111-01 

.414-02 


.259-00 

Ydr 

.299-03 

.152-02 

.656-04 

.618-05 

.245-01 


-.551-03 

L<Sr 

.342-01 

-.764-00 

.977-02 


-.209-01 


-.227-01 

N<Sr 

-.156-00 

-.172-00 

-.227-00 

-.127-01 



-.338+01 


Random Bias, 

Misalignments , and Center-of-Gravity Position Errors 

Parameter 


Yp 

Yr 

- y -p 

Ti 

e cgx 

. e cgz 

yb 

-.957-05 


-.278-03 


-.239-02 


.863-02 

lB 

-.109-01 

-.979-02 

-.368-01 



.397-00 

.287-01 

nB 

-.546-03 


-.323-02 

-.406-02 



-.179-02 

Lp 

-.658-03 


-.103-02 

.156-02 

.686-04 


.669-02 

Np 


-.195-04 

-.152-03 

.161-04 




Lr 

-.894-03 

-.188-02 

-.261-02 



. 109+00 


Nr 

-.162-03 


-.406-03 


-.150-02 

-.396-01 

-.994-02 

Y6a 

-.171-04 

.487-05 

.103-03 

-.709-06 

.772-03 

.138-01 


L<5a 



.135-01 

-.973-02 

.859-01 



N6a 



-.193-03 





YSr 

-.512-06 

- . 466-05 

.485-05 


.115-03 

.394-02 

-.464-02 

Lor 

-.351-02 

.174-02 

-.143-02 

.663-01 

.176-01 

-.616-01 

-.246-02 

NSr 

-.255-03 

.410-06 



.768-02 

.250-03 

.651-01 


-.439-03 


Mean Errors — Accelerometer Positions 

and Angl 

e-of-Sideslip Vane Position 

Parameter 

e ax 

e az 

e vx 

YB 

-.862-02 

-.863-02 

.158-05 

lB 

.397-00 

-.287-01 

.297-04 

NB 

-.166-01 

.179-02 

.614-04 

LP 

.117+00 

-.669-02 

-.156-04 

Np 

-.261-01 

.127-02 

.722-05 

Lr 

. 109+00 

-.321-01 

-.860-04 

Nr 

-.396-01 

.994-02 

.175-05 

Y<$a 

.138-01 

-.736-03 

-.131-05 

L6 a 

.510-01 

.236-02 

-.863-04 

N6a 

.546-01 

-.364-02 

.878-06 

Y<5a 

.394-02 

.464-02 

.368-06 

L6a 

-.616-01 

.246-02 

-.269-04 

N6a 

-.115-02 

.439-03 

-.573-04 
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Table 14 


DOMINANT RANDOM ERROR SOURCES FOR LATERAL PARAMETER IDENTIFICATION 


Parameter 

Biases Not 
Estimated 


Yg 

b ny’- e cgx* e cgz 

Y r * £ cgx * £ cgz 

YSa 

e cgx 

Y* , e , e 

r cgx* cgz 

Y<5r 

b ny ’ £ cgx * £ cgz 

£ cgx» £ cgz 

Lg 

bp ) Yr » e cgx 

Y r ’ e cgx 

Lp 

£ cgx 

£ cgx 

Lr 

Y., e , e- 

p cgx eg z 

Y. , £ , e 

p ’ cgx* cgz 

L6a 

b p> e cgx 

e* » y. . e 
p* y r > cgx 

L6r 

b P ’ £ cgx 

e* , Y* . Y* « £• 
p Y p ’ Y r * cgx 

Ng 

bp j b^’ , Y£ > £ cgx 

^r * £ cgx 

Np 

e cgx 

e cgx 

Nr 

bp , Yj. * £ cgx » e c.gz 

£ 9 £ 
cgx ' cgz 

N6a 

e cg X 

£ 

cgx 

N6r 

e r ’ Y r * £ cgx 

e. , y. 
r r 
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Table 15 


EFFECT OF INSTRUMENT REMOVAL ON PARAMETER DEVIATIONS FOR BASELINE SET 
OF ERRORS AND LATERAL EQUATIONS OF MOTION 


Parameter 


YB 

Y6a 

YSr 

LB 

Lp 

Lr 

L<5a 

L<$r 

NB 

Np 

Nr 

N6a 

N6r 


** 

Nominal 

All 

Instruments 

Value 


wmsm\ 

-.157s" 1 

.683-2 

.636-2 

-.00338s" 1 

.554-2 

.694-2 

•0246s 1 

.295-2 

.305-2 

-15.98s" 2 

.247 

.271 

-1.608s" 1 

.375-1 

.585-1 

.384s" 1 

.389-1 

.595-1 

10.92s" 2 

.941-1 

.764-1 

2.54s" 2 

.862-1 

.557-1 

6.563s" 2 

.232-1 

.206-1 

-.0997s" 1 

.490-2 

.132-1 

-.343s' 1 

.819-2 

.206-1 

.707s' 2 

.833-2 

.280-1 

-3. 902s~ 2 

.497-1 

.434-1 


No Angle of 
Sideslip Vane 

(B) 


Biases 


.692-2 .647-2 



Biases 


.720-2 


No Yaw Angle 
Accelerometer 

(r) 


Biases 


,654-2 .792-2 


585-1 .380-1 .597-1 .495-1 .671-1 


,553-2 .700-2 .548-2 .706-2 .553-2 

.295-2 .304-2 .283-2 .295-2 .292-2 

.250 .270 .270 .264 .339 

,380-1 .597-1 .495-1 .671-1 .407-1 

.407-1 .609-1 .143 .105 .418-1 

,967-1 .786-1 .204 .148 .101 

.887-1 .561-1 .227 .126 .830-1 

.240-1 .214-1 .367-1 .249-1 .250-1 

,491-2 .134-1 .510-2 .146-1 .476-2 

.864-2 .210-1 .896-2 .211-1 .899-2 

,864-2 .282-1 .168-1 .242-1 .125-1 

.520-1 .455-1 .513-1 .431-1 .905-1 


.899-2 

.125-1 

.905-1 


No B, P, 
or r 


Biases 


.694-2 
.803-2 
.298-2 
.416 

657-1 i .557-1 I .739-1 


.915-1 


358-1 .243-1 



Biases" means that all biases are error sources; "b<j> means all biases except b A are directly estimated. 


Dimension is : s = seconds 




















1 . 


Removal of the angle-of-gideslip vane has very little effect on the 
accuracy of the parameters. 

2. Removal of the roll angular accelerometer has a sizeable effect on 
the accuracy of Lp, Lr, L<5a, L6r, N$, and N6a. The largest effect 
is on Lr with the standard deviation increasing from about 10% to 
37% of the nominal parameter value. 

3. Removal of the yaw angular accelerometer has a marked effect on the 
accuracies of L3, N6a, and N6r, with Ndr's standard deviation increasing 
by 82%. 

4. Removal of the sideslip vane and the roll and yaw angular accelerometers 
simultaneously caused standard deviation of Y6a, Y6r, and Np to de- 
crease. The maximum increase was for Lr. 

5. Estimating biases caused little improvement in the identification accuracy 
More standard deviation magnitudes were increased than decreased, 

which again points out the fact that biases are not the prime source 
of error for lateral parameters. 
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IV 


SUMMARY AND CONCLUSIONS 

Two techniques - ensemble analysis and simulated data analysis - have 
been developed for determining the effects of instrumentation errors on 
identified aircraft parameter accuracy. The instrument errors contaminate 
the measured data which is used in the identification process. The analysis 
techniques are based on the assumption that an output error identification 
technique such as the modified Newton - Raphson algorithm is used. 

The two statistical analysis techniques have been coded into a digital 
computer program which allows rapid assessment of the error effects. The 
uses which can be made of this program include the following: 

1. The determination of the effect of instrumentation errors on the 
statistical accuracy of the stability and control derivatives and 
other parameters identified from flight test data can be made. This 
includes the mean error and standard deviation of each of the para- 
meters identified. The contribution of each error source on each 
parameter is determined. 

2. The effects of such variables as aircraft type and flight condition, 
control input sequence, and data sampling rate on the accuracy of 
the identified parameters can be determined. 

3. Tradeoff studies can be made between instrument quality and identi- 
fication accuracy. 

4. Different combinations of instruments can be studied for use in 
collecting the flight data. 

5. Tradeoff studies between fewer instruments with greater quality and 
more instruments with larger errors can be made. 

6. The necessary instrument accuracy required in a flight test program 

to allow identifying aircraft parameters to a desired level of certainty 
can be specified. 
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To exercise this program and determine the effect of some of the instru- 
ment errors on the variation of the identified stability and control deriva- 
tives, a study was conducted using the McDonnell F-4C aircraft. Both longitu- 
dinal and lateral equations of motion were utilized. Error sources assumed 
present were output measurement biases, scale factor errors, and correlated 
errors due to instrument misalignment, etc. Some of the general conclusions 
which can be made from the results of this study include: 

1. The instrumentation errors which are not due to noise can have 
very large effects on the identification accuracy. In the tests 
run, a growth in the standard deviations of many of the identified 
parameters exceeded an order of magnitude. 

2. Improvement can be made in the accuracy of the parameters by 
identifying the values of the dominant bias errors affecting the data. 
For the longitudinal equations of the F-4, identifying measurement 
biases generally improved the parameter accuracy obtained. However, 
identifying biases for the lateral equations was of little value 
because the dominant errors were from other sources. Adding the 
capability of identifying error sources which are of minor importance 
generally reduces the overall accuracy of the identified parameters 
because the same information is used to determine more quantities. 

3. The ratio of elements in the weighting matrix R of the cost function 
used by the Newton - Raphson scheme has an important roll in the over- 
all accuracy of the identified parameters. Changing the ratio of 
these elements changes the sensitivity of the identified parameter 
variations due to each error source. The effect of an error source 
which has a large variation can be minimized with the proper choice 

of the weighting matrix elements. 

4. The key output of the analysis program is the sensitivity matrix of 
the stability and control derivative errors to each of the instrument- 
ation error sources. With this matrix and an estimate of the instru- 
ment quality available, the test engineer can determine what the 
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accuracy of the identified aircraft parameters will be. He can 
also specify what instrumentation quality is necessary in order 
to identify the parameters to a desired level of accuracy. 

5. The control input sequence is a very important part of minimizing 
the effect of instrumentation errors and the uncertainty in the 
accuracy of the identified parameters. Also important to the 
parameter accuracy is the data sampling rate and the length of 
the data span. 

The above conclusions were based on linear analysis. Further conclusions 
can be expected if dynamic errors, control measurement errors, and known non- 
linear errors would also be included. 

Finally, it must be emphasized that no extensive laboratory or inflight 
study has been made of the kinds of errors that are prevalent in most flight 
instrumentation, including a statistical description of the error magnitudes. 
The range of instrument errors used in the study of this report were based 
on conversations with flight test personnel and reviewing instrument company 
literature. The information obtained from these sources was not in the statis- 
tical form necessary for the analysis* so assumptions had to be made. 
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V 


RECOMMENDATIONS 


In order to make general conclusions concerning the effect of instrument- 
ation errors on aircraft parameter identification, the following extensions to 
the study presented in Section III should be made: 

1. The effects of all known instrument errors not studied in this report 
should be determined. This includes instrument dynamics, control measure' 
ment errors, and non-linear error sources. 

2. Different types of aircraft and different flight regimes should be 
studied to determine their effect on the overall conclusions. 

3. A method of finding the "best" control input sequence corresponding 
to the available instruments used to collect flight data should be 
found. Always using the best control sequence in each study would 
remove the dependence of the accuracies on the control input. If 
the same control input is used when studying different sets of 
instruments, it is conceivable that one Instrument set will produce 
better results solely because the input sequence is more favorable 
to that set. 

4. Similarly, a method should be obtained for finding the most favorable 
weighting matrix used in the cost function of the Newton - Raphson 
identification algorithm. This matrix should tend to minimize the 
effect of the most prominant unidentified error sources. 

5. It should be established whether other error sources in addition to 
output biases can be identified. In cases where these other error 
sources are more dominant, identification accuracy can possibly be 
improved by identifying them as parameters and removing their effect. 

6. In order that instrumentation quality can be specified to meet the 
flight test objectives, a method must be established to define what 
constitutes an acceptable level of aircraft parameter uncertainty. 


77 



Based on the results of studying the F-4C aircraft, the following re- 
commendations are offered to insure that adequate instrumentation is provided 
for identification purposes: 

1. It was shown that instrumentation errors can produce significant 
increases in the identified parameter uncertainty. Therefore, 
the instrumentation quality required, including basic instrument 
accuracies, mounting accuracies, and recording equipment accuracies, 
should be specified prior to any flight test program. 

2. To insure that these specifications are met, laboratory and flight 
test studies should be made of all aircraft instrumentation used 
to collect data for identification. Statistical measurements of 
the instrument accuracy should be obtained to provide compatibility 
with the analysis methods. 

3. In addition to using instruments of necessary quality, care should 
be taken to align the instruments' sensitive axes, to calibrate the 
instruments, and to measure the instruments' and c.g. locations within 
the tolerances specified. 

This study has also shown that the specification of flight instrumentation 
may be sensitive to the particular flight control input sequence. Therefore, 
it is recommended that a means be developed to display the desired input to 
the pilot or to generate this input directly into the control system. Also, 
a method should be provided to tell the pilot if sufficient information is 
contained in the data collected during each maneuver which will enable identi- 
fication of the parameters to the level of accuracy desired. 


78 



APPENDIX A 

SUMMARY OF EQUATIONS CODED IN ANALYSIS PROGRAM 


A.l Explanation of Time Increments Used by Program 

The three time steps which are used in the program are: 

At - This is the time step of the Runge-Kutta integration 
L 

package. It is set small enough so the effect of the 

highest frequency dynamics is correctly simulated. 

This term is usually governed by the control or output 

measurement lag with the smallest time constant. 

At - This is the sample time increment used by the identifi- 
s 

cation process. When measurement lags are ignored, it 

is also used as the integration step. 

T - This is the time delay that governs when sampling of the 

control input is made after time zero. Normally, 

0 < T < At . Use of T enables sampling the control input 
— — s 

at different time points than the output measurements. 

A. 2 Simulated Data Analysis Subprogram - Mode 1 

A. 2.1 Simulated Airplane Equations 

The equations governing the simulated airplane are: 

x = Fx + Gu input ; x(0) = x Q (A.l) 


y_ = Hx + Du. 

'T input 

(true output) 


(A. 2) 
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y z - Ty T + B 

(indicated output - no lags) 

K ‘ *Vl + Vi ; V” ■ y i (o) 

Ceffect of lags) 


(A. 3) 


(A. 4) 


y i = 


Li 


+ w. 


(recorded output. If lags are present, this 

samples y every At seconds starting at 

time At s 

s 


(A. 5) 


A. 2. 2 Control Input Equations 


The equations governing the recorded control input are: 


U input 

(true control which is input to program) 


(A. 6) 


u = T u + B (A. 7) 

I c c 


(indicated control with no lags) 


u 


L 


u 


mi 


- f c u l + F C V u L (0) = Ul (0) 

(effect of lags) 

* 


U Li + 


w 


ci 


(A. 8) 


(A. 9) 


(recorded input. If lags are present, this samples 
u^ every At seconds starting at time T . 


* Note The control measurement noise is only used in the Monte Carlo 
Mode (Mode 2) of the Simulated Data Analysis. It is set 
to zero otherwise (Mode 1) . 
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A. 2. 3 Ensemble Analysis Subprogram Equations 


The random errors are first analyzed. To do this, the follow- 
ing equations are first integrated and evaluated in the ensemble 
analysis subprogram: 


x - Fx + Gu ; x(0) = x^ 

(At integration step) 



d_ 

dt 




( 0 ) 


1 


(for each initial condition estimated) 





/V 

X 


+ 



u 


(from Eq. A. 2) 


(A. 10) 


(A. 11) 


(A. 12) 


(A. 13) 




(for each initial condition estimated) 


(A. 14) 


3£ _ 

9 *b ' 


1 


(for each bias term estimated) 


W - 




9 ?il 


(A. 15) 


(A. 16) 
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n 


T 



I 




Invert 




(A. 17) 


(A. 18) 


A. 2. 4 Simulated Data Analysis Subprogram Equations for Random Errors 

Next, each of the random output measurement errors is 
sequenced. The object is to evaluate: 


aj 

3p 



i=l 


and 


6p 


8 2 J 

-1 

3J 

_ 3 P 2 . 


3p 


(A. 19) 


(A. 20) 


First, the effect of each unestimated output measurement bias 
is computed. This affects only the residual of Eq. (A. 19) as: 

( y i ~ = Y 1± ~ y i = B (A. 23) 

Then, each of the output transformation errors of T is sequenced. 
Again, the effect on Eq. (A. 19) is: 

( y i - y ± ) = T(Hx ± + Du i ) - (Hx i + Du ± ) = (T - I) (Hx. i + Du ± ) (A. 22) 
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Then, each of the unestimated initial condition errors is sequenced. 

This requires the integration of Eq. (A.l) with x(0) = Xq. Then, from Eq.(A.2), 


(y. - y.) = (y T . - y.) = h(x. - x.) 

i i *i i xi 


(A. 23) 


Next, the control measurement errors due to T and B are sequenced. 

c c 

To do this, the previous results of integrating Eq. (A. 10) are first 
stored. These state values represent the true airplane equations 
in this case. So 


x. = x. 
l i 


(A. 24) 


Then, Eqs. (A. 10), (A. 11), and (A. 12) are reintegrated using the At g in- 
tegration step and modified u .. Equations (A. 12) - (A. 18) are then 

mi 

reevaluated. For bias errors 

u . = u + B (A. 25) 

mi i c 


For scale factor errors 

u = T u (A. 26) 

mi c i 


For the random errors, the 

r 

E {6p6p T } - E noise + ^ 

J-l 

where r is the number of random error sources . E is the co- 

noise 

variance of the parameter estimate errors due to the output noise w^. 
Enoise e; *- t ^ er set to E< 1 * (A. 18) or computed using the Monte Carlo 
option with 


total parameter covariance is 

(5p.6p.^) (A. 27) 

3 J 



(A. 28) 


at each time point . 
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A. 2. 5 Mean Errors 


The mean errors are also sequenced. These include some 

elements in the T matrix (Eq. (A. 3) and all elements in F and F 

m c 

(Eqs. (A. 4) and (A. 8)). If there are no lag errors input, the evalua- 
tion is just like Eq. (A. 22). 


If there are control measurement lag errors, Eq . (A. 8) needs 

to be evaluated with u = u, over each interval At . The Integra- 

J- s 

tion step is At^. The stored is sampled every At g but begin- 
ning T seconds after time zero. 


If there are output lags, Eq. (A. 4) needs evaluation. This 
requires reintegration of Eq. (A.l), with x(0) = x Q , and Eq. (A. 4) of 
the form 


y T = y T + F „ ( Hx + Du ) 
L m L ci 


(A. 29) 


The step size is At . The output of Eq. (A. 29) is sampled every 

i-i 

At seconds with 
s 



(A. 30) 


Again, Eqs* (A. 10) - (A. 18) need to 
lags. Equations (A. 10) - (A. 12) use 
however. The effect on Eq. (A. 19) is 

‘ A.1 - <H *i + Du mi> 


be reevaluated in case of control 

an Integration step size of At , 

s 


(A. 31) 


A. 3 Simulated Data Analysis Subprogram - Mode 2 


The random errors contained in B, part of T, B T , and x are 
generated at the beginning of each run using the input standard devia- 
tions and the random number generator. They are held constant during 
each single Monte Carlo run. The random noise w^ and w c ^ are regener- 
ated each sample point during each run. Each of the mean errors in T 
plus elements of F^ and F^ are set equal to the values input and are 
not changed during any of the runs. 
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For output measurement errors in T, B and w^, only, the 
residual (y^ - y^) in Eq. (A. 19) is computed by 

y. - y 4 = T (Hx + Du.) + B + w. - (Hx. + Du.) (A. 32) 

■'i i i i 11 i 

For random initial conditions added, Eq. (A.l) must be integrated 
each time and Eq. (A. 32) gets changed to 

y i - y = T(Hx i + Du ± ) + B + w ± - (Hx + D Uj ,) (A. 33) 

For non-lag control measurement errors, Eqs. (A. 7) and (A. 9) 
get combined so that at each sample point 


u,=Tu + B + w . 
mi c c ci 


(A. 34) 


where is randomly generated each time point. Because of this 
change, Eqs. (A. 10) - (A. 12) require integration each pass through, 
and Eqs. (A. 13) - (A. 18) require re-evaluation each pass through. With 
these changes, Eq. (A. 33) becomes 


y - y ± = T(Hx i + Du ± ) + B + w ± - (Hx ± + Du^) 


(A. 35) 


For lags in the control input measurements and output measure- 
ments, two integration step sizes have to be used. First, Eqs. (A.l) - 
(A. 4) are evaluated using the At T step. The resulting y is sampled 

Li L 

every At seconds starting at time At . Then, Eq. (A. 5) is evaluated 
s s 

and the resulting y^ stored. Next, Eqs. (A. 7) - (A. 8) are evaluated using 

the At step. The resulting u is sampled every At seconds starting 
L Lj S 

at time T. Then, Eq. (A. 9) is evaluated and the resulting u . stored. 

mi 

Next, Eqs. (A. 10)- (A. 12) are evaluated using the At step and the stored 

u .. Then, Eqs. (A. 13) - (A. 18) are evaluated. Then, the residual for 
mi 

Eq. (A. 19) is computed using 


A 



(Hx + Du.) 

mi 


(A. 36) 
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In Mode 2, w . is used lust like Wj . 

ci i 

In Mode 2 , several Monte Carlo runs are made. In 

each run, the values of the random elements in T, B, and x are ran- 

o 

domly generated at the beginning of the run. For comparison, the 
standard deviations used to generate these values are the same as 
those values used one at a time in Mode 1. These random values 
remain constant throughout a run, but change from run to run. The 
white noise vector w^ is randomly generated, and it changes every 
sample point throughout the run. Those mean error sources in T, 

F m’ and F c which aren 't random are held constant throughout all the 
runs. The error Ap.. in the parameter vector obtained from each run 
is saved. For m Monte Carlo runs, the mean error in p is 


Ap A E {Ap} - ^ J Ap 

j=l 


(A. 37) 


The standard deviation about this mean is 

m 

E {6p6p T } = ~ J (Ap - Ap) (Ap - Ap) T (A. 38) 

3=1 


For small values of the error sources, the results of Eq. (A. 38) should 
match those of Eq. (A. 27). Also, Eq. (A. 37) should match the value of 
the mean obtained in Mode 1. 
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A. 4 Ensemble Error Analysis Subprogram 


A. 4.1 Output Errors 

The basic equation which is used to compute changes in the 
parameters is 


6p = - 


r 3 2 J~ l ~ 1 

L 3p 2 J 


3J 

3p 


= + 


n / T _i n 

v K R -i Ox' 

) 3p 3p 

L i-1 ( 


2 ** 


3p 


--1 f ^ \T 

R (y i _y i ) 


i=l 


(A. 20) 


The are the sampled measurements taken from the aircraft 
(Eqs. (A.l) - (A. 5)) and the £ are the simulated values obtained from 

the measured control input (Eqs. (A. 6) - (A. 9) and Eq. (A. 10)). The 
term 3y^ comes from Eqs. (A, 11) - (A. 16). 

3p 

It is assumed that if lag errors exist, lag terms are present 
in all control input and system output measurement equations. If the 
measurements of y^ have lags, the dynamic equations can be written as 

x = Fx + Gu (A.l) 

y T = -F y T + F (T [H x + D u] + B ) (A. 39) 

L m L m 


Otherwise, the output equations are of the form 
y = T(Hx + Du) + B 


(A. 40) 
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biases, 


The sensitivity of y_^ to each error e is computed. For 


h. = i 
3e 


For errors due to the matrix T, the sensitivity can be found as 


»» 


For initial condition errors, 


9 £ 

3e 


3e 


where the term 


3x 

3e 


comes from integrating 


d_ 

dt 





I 


When output lags 
equation comes from Eq. 


are present, 
(A. 39), and 


the sensitivity differential 
is 


d_ 

dt 



y + Hx + Du 
Li 


F 



m 3e 


The results of Eqs. 
8y i 

- — for each error 
3e 


(A. 41) through (A. 45) are combined into a general 
e which affects the output measurements. 


(A. 41) 


(A. 42) 


(A. 43) 


(A. 44) 


(A. 45) 
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A. 4. 2 Control Input Errors 


With respect to Eq. (A. 20), the error effect on the term ^ 
is computed. The differential equations of the simulated aircraft are 

"l = _F c U L + F c (T c u + V 5 V°> = ^CO) 
x = Fx + Gu 


Li 


y . = Hx . + Du T . 
l l Li 


For control measurement lags the following equations, that 
from Eqs. (A. 46 and (A. 47), are first integrated. 


come 


d_ / 3x 
dt \ 3e 


a. (h, 

dt \ 3e J 


3u T 

= F lf +G JT > 


\ 3u 3F 

) = ~ F cjr + dT + 


3x 


u) 


(A. 46) 
(A. 47) 
(A. 48) 


(0) = 0 


(A. 49) 

3u l 

3u 

(A. 50) 

— (0) = 

3e 



Then, from Eq. (A. 48) 

3y . 3x. 3u 

X i T -i 

sT ■ H ar + D 5^ 

For control scale factor errors, the program first integrates 

8T^ 

11 • 

de 

Then, the sensitivity is computed as 


, # (0) = 0 


3y 34. 

_i = H 1 


+ n 9T 
+ D c u. 

3e 1 


3e 3 e 

For control biases, the program integrates 


(A. 51) 


(A. 52) 


(A. 53) 
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9_ 

at 



F 


3x 3B 

— + G — c_ . 
ae 3e * 



0 


(A. 54) 


Then, the sensitivity becomes 





(A. 55) 


A. 4. 3 Effect on Parameter Estimates 


The sensitivity of parameter estimates due to output measure- 
ment is of the form 




3 2 J 

3p 2 


-1 


I 


i=l 


3y. . 3y . 

_JL R -1 
3p * 3e 


(A. 56) 


3 y< 


where comes from evaluating Eqs. (A. 41) through (A. 45). 


The sensitivity of parameter estimates due to control input 
errors is of the form 


3e 


(<S P ) 


3 2 J 


-,-1 n 


L^p 2 J 


i=l 




where 


3e 


comes from evaluating Eqs. (A. 51), (A. 53), and (A. 55). 


(A. 57) 


For random errors, the total parameter covariance is 

E (dp dp T ) = E no . se + ± (6p)E {e e T } J- (6 p ) T 

For mean errors, the mean parameter error is 
E (dp) = — ^P-^£{e} 


(A. 58) 


(A. 59) 
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